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ABSTRACT 


In  the  Stanford  Radio  Acoustic  Sounding  System  (RASS),  an  electro- 
magnetic signal  is  made  to  scatter  from  a moving  acoustic  pulse  train. 
Under  a Bragg-scatter  condition,  maximum  electromagnetic  scattering 
occurs.  The  scattered  radio  signal  contains  temperature  and  wind 
information  as  a function  of  the  acoustic  pulse  position. 

In  the  theoretical  work  on  RASS  to  date,  the  effects  of  such 
atmospheric  parameters  as  turbulence,  humidity,  mean  temperature,  and 
mean  wind  fields  on  the  propagating  acoustic  pulse  train  have  been 
ignored.  By  neglecting  these  parameters,  the  quantitative  analyses 
have  assumed  that  the  acoustic  wavefronts  act  as  large  perfect  spherical 
reflectors.  In  this  investigation,  RASS  performance  is  assessed  in  a 
real  atmosphere  where  "coherency"  of  the  acoustic  pulse  is  degraded  as 
it  propagates  vertically  into  the  lower  atmosphere.  The  only  assumption 
made  is  that  the  electromagnetic  wave  is  not  affected  by  stochastic 
perturbations  in  the  atmosphere. 

Coherency  of  vertical  acoustic-wave  propagation  is  described 
through  a perturbation- theoretic  method  and  Feynman’s  diagrammatic 
technique.  One  of  the  most  important  attributes  of  this  analysis  is 
that  it  systematically  and  explicitly  accounts  for  multiple  scattering 
of  acoustic  waves  in  the  presence  of  atmospheric  fluctuations.  The 
coherency  results  are  then  used  to  evaluate  the  strength  of  the  scat- 
tered electromagnetic  signal  from  the  acoustic  pulse  train  while  taking 
into  account  the  presence  of  turbulence,  mean  temperature  gradients, 
and  mean  wind  fields. 

It  is  concluded  that,  for  acoustic  pulses  with  carrier  frequencies 
below  a few  kilohertz  propagating  under  typical  atmospheric  conditions. 
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turbulence  has  little  effect  on  the  strength  of  the  received  radio 

signal  at  heights  up  to  a few  kilometers.  This  result  implies  that 

focusing  of  RF  energy  by  the  acoustic  wavefronts  is  primarily  a function 

_2 

of  sound  intensity  which  decreases  as  x , where  x denotes  altitude. 

The  effect  of  mean  vertical  wind  and  mean  temperature  on  the 
strength  of  the  received  signal  is  also  demonstrated  to  be  insignifi- 
cant. Mean  horizontal  winds,  however,  shift  the  focus  of  the  reflected 
electromagnetic  energy  from  its  origin,  resulting  in  a decrease  in 
received  signal  level  when  a monostatic  RF  system  is  used.  For  a 
bistatic  radar  configuration  with  space-diversified  receiving  antennas, 
the  shifting  of  the  acoustic  pulse  makes  possible  the  remote  measure- 
ment of  the  horizontal  wind  component. 
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Chapter  I 
INTRODUCTION 

During  the  past  few  years,  a ^adio-Acoustic-^ounding-System  (RASS) 
has  been  developed  for  real-time  temperature  profiling  the  first  few 
kilometers  of  the  lower  troposphere  [Marshall,  1972].  The  basic  physics 
of  the  RASS  is  as  follows. 

The  electromagnetic  refractive  index  of  air  (above  30  MHz)  in  the 
lower  troposphere  is  not  exactly  unity  but  is  a function  of  pressure, 
temperature,  and  humidity  [Bean  and  Dutton,  1968].  Because  of  its 
pressure  dependence,  this  refractive  index  can  be  altered  by  a short 
pulse  of  sound  from  an  acoustic  source.  An  electromagnetic  RF  signal 
generated  by  a radar  passes  through  the  sound  pulse  and  is  scattered  as 
a result  of  the  induced  refractive-index  variations.  The  scattered  RF 
signal  is  collected  by  a receiver  and  is  processed  to  determine  atmo- 
spheric temperature  as  a function  of  the  acoustic-pulse  position  (range). 
In  a static  atmosphere,  the  received  signal  is  maximized  when  a Bragg 
condition  is  established  between  the  acoustic  and  electromagnetic  sig- 
nals. This  occurs  if  the  electromagnetic  wavelength  is  twice  the 
acoustic  wavelength  and  results  in  an  in-phase  addition  at  the  receiver 
of  the  electromagnetic  signal  scattered  from  successive  acoustic  wave- 
fronts  in  the  pulse. 

A doppler  radar  measures  the  speed  of  the  sound  pulse.  Because  the 
speed  of  sound  in  air  depends  on  its  temperature,  the  temperature  pro- 
file of  the  lower  atmosphere  can  be  obtained  with  the  RASS  [Frankel  and 
Peterson,  1976].  Fig.  1 illustrates  the  geometry  of  a bistatic  RASS 
wherein  the  acoustic  source  is  located  between  separated  RF  transmit  and 
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Fig.  1.  REFLECTION  OF  ELECTROMAGNETIC 
ENERGY  FROM  AN  ACOUSTIC  PULSE  FOR  A 
RASS  B I STATIC  GEOMETRY. 

receive  antennas.  In  a monostatic  RASS,  the  antennas  and  acoustic  source 
are  located  at  the  same  point. 

In  an  early  study  at  Stanford  University  , the  reflecLiori  of  electro 
magnetic  energy  from  an  acoustic  pulse  propagating  vertically  into  the 
lower  atmosphere  was  analyzed  at  a carrier  frequency  of  85  Hz  [Marshall, 
1972].  In  that  study  Marshall  assumed  a static  atmosphere  ( no  turbu- 
lence and  winds  ),  and  he  made  no  attempt  to  determine  the  feasibility 
of  operating  the  RASS  at  acoustic  frequencies  higher  than  85  Hz  which 
offers  the  advantage  of  a smaller  system  and  an  increase  in  the  resolu- 
tion of  the  measured  atmospheric  parameters. 

In  contrast  to  Marshall's  analysis,  this  investigation  considers 
the  interaction  of  electromagnetic  and  acoustic  waves  when  the  RASS 
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is  opera ,ed  in  a stochastic  environment  characterized  by  turbulence, 
winds,  and  mean-temperature  gradients.  These  atmospheric  parameters  are 
important  when  evaluating  RASS  performance  because  they  could  affect 
the  spatial  coherency  in  and  between  the  acoustic  wavefronts  of  the 
transmitted  acoustic  pulse.  This  decrease  in  coherency  would  reduce  the 
received  signal  levels. 

An  analysis  of  the  effects  of  turbulence  on  the  coherency  of  the 
acoustic  wave  shows  that  the  acoustic  wave  amplitude  at  an  altitude  x is 
damped  by  a factor  exp[- <(^| {x,y^)  | ^ ^/2]  because  of  complex-phase 
intereference.  In  this  damping  factor,  4'i(x,yj^ ) is  the  first-order  pertur- 
bation in  the  complex  phase  of  the  acoustic  wave  due  to  turbulence.  As 
shown  herein  , these  complex-phase  perturbations  of  the  acoustic  wave 
are  small  in  the  lower  troposphere;  consequently  , for  acoustic  pulses 
with  carrier  frequencies  below  a few  kilohertz  propagating  (under  typical 
atmospheric  conditions)  to  altitudes  of  a few  kilometers,  turbulence 
has  little  effect  on  the  strength  of  the  received  radio  signal.  More 
particularly,  it  is  shown  in  Chapter  IV,  that  for  an  altitude  x of  the 
acoustic  pulse  much  less  than  the  coherence  length  of  propagation  x^, 
the  effect  of  turbulence  on  RASS  performance  is  insignificant.  This 
coherence  length  of  propagation  derived  in  this  investigation  is  given  by 


X = 

c 


54.8  (1.833  - m^) 

"c  ' k'  a'/3^ 

no  a 0 


6/(ll-6  m^) 


(1.1) 


In  the  above  equation  m^  is  a factor  determining  the  acoustic 

refractive-index  structure  parameter  in  the  presence  of  inhomogeneous 
2 

turbulence,  is  the  refractive-index  structure  parameter  at  altitude 
^o’  *^a  acoustic  wave  number,  and  A is  a proportionality  constant 


3 


( 

( 

i 


' measuring  the  outer  scale  of  turbulence.  It  is  observed  that  the  coheren- 

ce length  varies  inversely  with  frequency  and  the  strength  of  the  turbu- 
lence. As  an  example  of  the  application  of  Eq.  (1 . 1 ), consider  a RASS 

f 

operating  at  an  acoustic  frequency  of  f^  = 2 kHz  and  assuming  the  folio- 

2 “6  -2/3 

. wing  typical  values  for  the  parameters  in  the  equation  ; = 10  m , 

X = 1 m,  A = 2 , m = 1.33, and  the  velocity  of  sound  c = 340  m/sec, 

^ 0 Q 0 

I the  coherence  length  is  found  to  be  x^  = 4 x 10^  m.  This  value  of  cohe- 

rence length  implies  that  the  coherency  of  an  acoustic  wave  is  not 
affected  by  turbulence  when  f,=  2 kHz  and  x«  4 x lO^m. 

a 

[ Quiescent  atmospheric  conditions  such  as  temperature  gradients 

I cause  the  dispersion  of  acoustic  wavetrain.  Under  these  circumstances,  a 

I match  between  electromagnetic  and  acoustic  waves  can  be  obtained  by  the 

I use  of  a modified  Bragg-scatter  condition.  This  match  depends  on  the 

I temperature  difference  over  the  length  of  the  acoustic  wavetrain.  This 

I study  shows,  however,  that  typical  temperature  gradients  in  the  lower 

atmosphere,  which  are  on  the  order  of  - 6.5°K/km  , have  a negligible 
effect  on  the  received  power  of  a RASS.  The  effect  of  vertical  winds 
on  the  strength  of  the  received  radio  signal  can  be  overcome 
similarly. 

This  investigation  also  shows  that  mean  horizontal  winds  shift  the 
focus  of  reflected  electromagnetic  energy  from  its  origin,  resulting  in  a 
decrease  in  received  signal  level  when  a monostatic  RF  system  is  used. 
Because  the  acoustic  wavefronts  act  as  large  spherical  reflectors,  the 
principle  of  specular  reflection  , however,  can  be  utilized  in  a bistatic 
radar  geometry  to  remotely  measure  this  wind  component  ( in  addition  to 
atmospheric  temperature  profiles  ).In  this  RASS  configuration,  the  radar 
antennas  and  acoustic  sources  are  aligned  in  the  wind  direction  --  the 

i 
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transmit  antenna  upwind  from  the  acoustic  source  and  the  receive  antenna 
downwind.  This  technique  offers  the  promise  and  potential  of  measuring 
horizontal  wind  velocities  of  magnitude  up  to  a few  tens  of  meters  in  the 
lower  troposphere.  Consistent  with  these  theoretical  predictions  Frankel 
et  al.,[1977]  have  utilized  this  principle  of  specular  reflection  to 
measure  horizontal  winds  in  the  lower  troposphere. 

The  analysis  in  this  study  was  restricted  to  the  lower  troposphere  ; 
however,  it  can  be  modified  to  study  this  diffraction  phenomenon  in  any 
media.  The  mathematical  and  physical  representations  of  the  acoustic 
refractive-index  field  and  atmospheric  turbulence,  following  Tatarskii 
[1961,  1971]  and  Yaglom  [1962],  are  presented  in  Chapter  II.  Acoustic- 
wave  propagation  in  a turbulent  medium  is  discussed  in  Chapter  III, 
and  Feynman's  diagrammatic  approach  is  used  in  Chapter  IV  to  derive 

1 

measures  for  the  coherency  of  a vertically  propagating  acoustic  wave.  j 

Results  obtained  in  Chapter  IV  are  applied  in  Chapter  V to  develop  1 

asymptotic  expressions  for  the  electromagnetic  energy  scattered  from  an  < 

acoustic  wavetrain  perturbed  by  atmospheric  turbulence,  mean  temperature, 
and  mean  wind-fields.  Chapter  VI  discusses  the  effect  of  the  mean  wind 
field  on  received  doppler  frequency.  Chapter  VII  summarizes  the  conclu- 
sions of  this  research. 
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Chapter  II 

ACOUSTIC  REFRACTIVE-INDEX  FIELD 


Acoustic  and  electromagnetic  waves  propagate  in  a stochastic  medium; 
however,  the  naturally  induced  fluctuations  in  the  acoustic  refractive 
index  are  approximately  one  thousand  times  greater  than  those  in  the 
electromagnetic  index  [Little,  1969].  For  this  reason,  natural  fluctua- 
tions in  the  electromagnetic  refractive  index  are  ignored  in  this  study. 

The  acoustic  refractive  index  depends  on  wind,  temperature,  and 
humidity,  and  the  stochastic  variations  in  these  parameters  constitute 
the  phenomenon  known  as  "atmospheric  turbulence."  In  this  chapter, 
random  fields  are  described;  the  mathematical  analysis  follows  Tatarskii 
[1961,  1971]  and  Yaglom  [1962].  This  analysis  ’s  then  extended  to 
describe  atmospheric  turbulence  and  its  effect  on  acoustic  refractive- 
index  fields  in  the  lower  troposphere. 

A.  Random  Fields 

Physical  quantities  such  as  atmospheric  temperature,  wind,  density, 
and  pressure  are  continuously  distributed  in  space  and  time  and,  conse- 
quently, defy  deterministic  description.  These  can  be  described,  how- 
ever, in  terms  of  random  fields  which  are  random  space-time  functions. 

In  this  section,  random-field  theory  is  described,  with  emphasis  on 
meteorological  applications. 


1 . Statistical  Representation  of  Homogeneous  Random  Fields 

Let  p(r)  be  a random  field,  where  r is  the  space  coordinate 
and  let  denote  the  mathematical  expectation  of  a random  function. 
The  mean  of  the  random  field  is  ^p(r)^,  and  its  covariance  is  defined 
as  , 


FMBim  iiw  Mat  numo 


Bp(r,r’) 


^ P(r)  - <^p(r))>]  [p*(r')  - <(pV')^]^  (2.1) 


(asterisks  denote  complex  conjugates).  A random  field  is  "homogeneous" 
if  its  probability  distribution  function  is  invariant  under  a space 
translation;  it  is  "homogeneous  in  the  wide  sense"  if 


and 


^P(>^)^  = constant 


Bp(r,r')  = Bp(r-r',0)  = Bp(r-r') 


(2.2a) 

(2.2b) 


Henceforth,  such  fields  will  be  referred  to  as  "homogeneous  random 
fields . " 

For  an  isotropic  homogeneous  random  field,  B (r)  = B (r) 

P P 

where  jr|  = r is  the  magnitude  of  the  vector  r.  When  <(^p(r)^ 
p(r)  has  a Fourier-Stiel tjes  expansion. 


= 0, 


such  that 


P(>")  =JJJ' e^'^''^  dW(K:) 


<^dW(K)^  = 0 

^dW(K)dW  (k')^  = 6(<-k:' )$  (K)d^K  d^' 


and  <l'p(K)  >0  is  the  spectral  density  of  the  field  where 


(2.3a) 


(2.3b) 


Bp(r)  -JJJ 'l>p(ic)d^ 


(2.4) 


When  the  medium  is  isotropic,  4>p(K:)  = ^p('^) 


1 


--.i 

1 


I 

I 

I 


The  above  equation  describes  globally  homogeneous  random 
fields;  however,  in  the  lower  troposphere,  meteorological  variables 
are  not  statistically  homogeneous.  Homogeneity  is  violated  by  large- 
scale  events  and,  therefore,  it  is  necessary  to  consider  locally  homo- 
geneous random  fields. 

2.  Locally  Homogeneous  Random  Fields  and  Structure  Function 

A random  field  is  locally  homogeneous  if  the  probability 
distribution  of  [p{r)  - p(r')]  is  invariant  under  a joint  space 
translation  of  the  points  7 and  7'.  Such  fields  are  described  by 
"structure  functions"  first  introduced  by  Kolmogorov  [1941].  This 
function,  for  a real-valued  random  field,  is  defined  as 

2 

Dp(7,7')  =<^j[p(r)  - <^p(r))']  - [p{7')  - <^p(7’))>]j  y {2.5a) 

Physically,  Dp(r,r')  depends  on  the  scale  sizes  of  the  inhomogeneities 
in  the  atmosphere,  which  are  less  than  |7  - 7' | . 

For  a locally  homogeneous  random  field, 

Dp(r,7')  = Dp (7-7’ ,0)  E Dp(7-7')  (2.5b) 

If  the  field  is  also  locally  isotropic,  then  Dp(7)  = Dp(>")- 

The  Fourier-Stieltjes  expansion  of  a locally  homogeneous 
random  field  is 

CO 

p(r)  = p(0)  + Ip-?  - 1]  dW(7)  (2.6) 

• CO 


where  p(0) 


is  a random  variable  and  ap  is  a random  vector. 


If 


J 


I 

I 


= 0 and 
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p(r)  = p(0)  - 1]  dW(K) 


(2.7a) 


such  that 


<CdW(K)>  = 0 

^dW(K)dW  (ic’)^  = 6(ic-ic' )'I'p(ic)  d^  ic  d^ 
'I>p(tc)  > 0 


(2.7b) 


(2.7c) 


(2.7d) 


Dp(r)  = 2 JJJ'd  - cos(tc- r)]'t'p(<)  d^< 


(2.8) 


If  this  field  is  also  isotropic,  then  ♦p(^)  = '*’p('^)-  7he  covariance 
and  structure  function  of  a homogeneous  and  isotropic  random  field  are 
related  by 


Dp(r)  = 2[Bp(0)  - Bp(r)l 


(2.9a) 


Because  Bp(“)  0 a physical  random  field,  Dp(®)  = 2Bp(0)  and 


B.(r)  - i [D  (-)  - D (r)] 


(2.9b) 


In  this  subsection,  a three-dimensional  spectral  expansion  of 
a locally  homogeneous  random  field  was  made,  and  a two-dimensional 
spectral  expansion  is  described  below. 


3.  Two-Dimensional  Fourier-Stiel tjes  Representation  of  Local  1' 


Homogeneous  Random  Fields 


The  two-dimensional  Fourier-Stiel tjes  expansion  of  a random 
field  p(rj  in  the  plane  x = constant  is 


10 


p(r)  = p(x,0,0)  + - l]  dSp(x,<^)  (2. 10a) 


where 


= (<2.1^3) 


(2.10b) 


- (y.z) 


(2.10c) 


p(x,0,0)  is  a random  function  and 


^ds  (x,K  )ds*(x',<')>  = 6(k^-ic' )F„(x-x' ,K  ) d^  < d^  k'  (2.10d) 
'PaP  a ''  a a p a a a 


) - F (-x,ic  ) 
Papa 


(2.10e) 


As  shown  by  Tatarskii  [1971],  the  two-dimensional  spectrum 


and  the  structure  function  are  related  by 


Dp(r-F-)  - Dp(x-x',0,0)  = 2jfy*[l  - cos  • (y^-y;)]  Fp(x-x' ,<^)d^ 


(2.10f) 


When  5 = x-x'  = 0,  y-y'  = n,  and  z-z'  = c,  the  above  equation 


becomes 


Dp(0,n. rj  - 2 - cos(K2ri+K3C)]  Fp(0,iC2,iC2)  ^<2  (2.10g) 


The  two-  and  three-dimensional  spectral  densities,  defined  in  Eqs. 


(2.10a),  (2.10d),  (2.10e),  and  (2.4),  respectively,  are  related  by 
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In  addition. 


Op(r)  . - Fp(x,.^)cos(K„  y_^)]  (2.12) 

-oo 

and 

OO 

—CO 

The  two-dimensional  Fourier-Stiel tjes  representation  of  locally  homo- 
geneous random  fields  will  be  used  in  Chapters  III  and  IV  to  describe 
the  coherency  of  a propagating  acoustic  wave. 

In  the  spectral  expansion  of  random  fields,  it  has  been 
assumed  that  the  medium  has  constant  mean  characteristics.  In  the 
lower  troposphere,  however,  the  intensity  of  fluctuations  is  smoothly 
varying,  and  the  above  descriptions  of  random  fields  must  be  modified. 


4.  Medium  with  Smoothly  Varying  Mean  Characteristics 

Following  Tatarskii  [1971],  if  the  random  field  p(r)  is 
smoothly  varying,  then 


The  function  Dp°^(r-r')  represents  the  local  variation  of  fluctuation 
intensity  and  is  defined  for  arguments  with  magnitudes  much  smaller  than 
the  largest  scale  size  events.  The  function  Cp{(r+r')/2)  represents 
the  smooth  variations  and  is  approximately  constant  when  |r-r'|  is  much 
smaller  than  the  largest  scale  size  event.  Tatarskii  demonstrated  that 


m 


P\ 


and  that 


(2.15) 


(2.16) 


These  equations  will  be  used  to  describe  a spherical  wave  traveling 
vertically  into  the  lower  atmosphere  which  is  locally  homogeneous  and 
has  smoothly  varying  mean  characteristics. 

B . Atmospheric  Turbulence 

The  physical  phenomenon  associated  with  the  deviations  of  energy 
transportation  processes  in  the  atmosphere  from  its  deterministic 
description  is  called  "atmospheric  turbulence."  The  medium-describing 
variables  are  random  functions  of  four  coordinates,  space  and  time. 

In  this  section,  the  effects  of  turbulence  on  the  acoustic  refractive 
index  are  considered.  Let  the  acoustic  refractive- index  n(r')  be 
defined  as 

n(r)  = <n(r))>  + (r)  (2.17) 


where  <\n(r)^  = 1 and  ^n^(r)^  = 0.  In  addition,  ni(r)  is  the 
fluctuating  component  of  n(r)  and  is  usually  much  smaller  than  unity. 


The  random  medium,  the  lower  part  of  the  troposphere,  generally  is 
not  statistically  homogeneous;  however,  the  fluctuating  acoustic  refrac- 
tive-index field  is  locally  homogeneous.  As  discussed  in  subsection  A. 2, 
such  fields  can  be  described  in  terms  of  the  structure  function  of  n(T')» 
defined  as 

Dn(r,r')  = <^{n■^{r)  - n^(7')}^)>  |7-r' | < (2.18) 


In  ■'■his  equation,  <(n(r)^  is  approximately  constant  for  |7-r' | < L^, 

where  is  called  the  "outer  scale  of  turbulence"  and  is  the  least 

distance  between  coordinates  r and  r'  for  which  B (7,r')  = 0. 

n 

For  points  separated  spatially  by  distances  greater  than  L^,  the 
covariance  function  is  zero. 

The  structure  function  for  locally  homogeneous  isotropic  turbulence 
was  first  obtained  by  Kolmogorov  [1941];  according  to  his  theory. 


Dn(r) 


n 


r <<  £ 

0 


a <<  r « L 
0 0 


(2.19) 


2 

where  is  the  "structure  parameter"  and  is  a measure  of  the  inten- 
sity of  fluctuations.  The  quantity  is  called  the  "inner-scale  of 

turbulence,"  and  the  interval  between  £ and  L is  the  inertial 

0 0 

subrange.  Energy  is  injected  through  some  outside  mechanism  into  L^. 
According  to  Kolmogorov's  cascade  theory  of  turbulence,  this  energy  is 
then  successively  transferred  to  smaller  scales  of  turbulence  until  it 
reaches  the  size  of  9.^  where  the  viscous  forces  dominate  and  energy 


I 
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is  converted  into  heat.  The  scales  of  turbulence  between  £ and  L 
^ 0 0 

I are  also  often  referred  to  as  "eddies"  [Tatarskii,  1971].  Generally, 

[ and  L„  are  functions  of  the  vertical  distance  x above  the 

[ 0 0 

ground.  According  to  Lawrence  and  Strohbehn  [1970],  L is  on  the 

i ° 

I order  of  one  to  hundreds  of  meters,  depending  on  the  altitude;  Gray  and 

Waterman  [1970]  show  to  be  on  the  order  of  a few  millimeters.  A 
model  for  the  outer  scale  of  turbulence  as  a function  of  altitude  x is 

i Lq(x)  = Ax^/^  (2.20a) 

I where  A varies  from  1 to  10,  with  a typical  value  of  2 [Taylor,  1968]. 

f Because  only  limited  measurements  exist  for  the  inner  scale  of  turbulence 

[Lawrence  and  Strohbehn,  1970]  and  for  mathematical  simplicity,  the  func- 
tional dependence  of  the  inner  scale  as  a function  of  altitude  is  modeled 
as  the  outer  scale.  Therefore, 

«.p(x)  = Bx^^^  (2.20b) 

where  A/B  = 10^  to  10^.  For  example,  at  x = 1000  m and  B = lO"^, 

the  size  of  the  inner  scale  of  turbulence  is  = 3 mm. 

0 

In  the  lower  troposphere,  the  refractive- index  field  is  locally 

homogeneous  and  has  smoothly  varying  mean  characteristics  in  the  vertical 

direction  x [Brown  and  Keeler,  1975].  The  refractive- index  structure 
2 

parameter  , as  defined  in  Eq.  (2.19),  is  a function  of  altitude  x. 

' Consequently, 

-m 

1.33>n,,  >0  (2.200 
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2 2 2 
where  ^ increases,  C^(x)  decreases  and  has  a 

maximum  value  near  the  ground.  According  to  Brown  and  Keeler  [1975], 

2 -6 

m^  is  approximately  equal  to  1.33, and  ranges  from  lO’  to 
10‘V^^^  at  Xq  = Im. 

C.  Refractive-Index  Spectrum  in  the  Lower  Troposphere 

For  isotropic  homogeneous  turbulence  in  the  inertial  subrange, 
Kolmogorov's  two-third  structure  function  law  results  in  the  spectrum, 

$^(ic)  = 0.033C^  p- « k « (2.21a) 

^0  0 

For  wave  numbers  k >>  21^/*.^,  is  small  and  the  corresponding 

range  of  the  structure  function  given  in  Eq.  (2.19)  for  r << 

Tatarskii  [1971]  used  the  following  form  for  the  spectrum. 

4)^(k)  = 0.033C^  e ^ ^ ^ (2.21b) 

0 

where  k = 5.92/J,  = 2^^/^  . 
m 0 0 


A disadvantage  of  the  above  representation  for  is  that  it  has  a 

singularity  at  k = 0.  As  a result  Tatarskii's  spectrum  does  not 
possess  an  autocovariance  function.  An  alternative  representation  of 
the  spectrum,  as  suggested  by  Strohbehn  [1968]  and  Lutomirski  and  Yura 
[1971],  is 


2 

0.033C  e 


" 7 2 ‘/Y-'2J176 


(2.22) 


where  k = 2v-r  and  L = 2ift  . 
0 0 0 0 


This  spectrum  is  flat  for 


ii 
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In  the  lower  troposphere,  the  turbulent  medium  generally  has  a 
smoothly  varying  mean  characteristic;  therefore,  substituting  Eqs. 


in  the  plane,  x = constant.  This  spectral  expansion  of  the  acoustic 
refractive-index  field  will  be  used  in  subsequent  chapters  to  describe 
acoustic-wave  propagation  in  a turbulent  atmosphere. 
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Chapter  III 

ACOUSTIC-WAVE  PROPAGATION  IN  THE  LOWER  ATMOSPHERE 

In  this  chapter,  the  basic  equations  for  acoustic-wave  propagation 
are  derived.  The  theory  developed  by  Tatarskii  [1961,  1971]  is  valid 
for  acoustic  wavelengths  much  smaller  than  the  inner  scale  of 

a 

turbulence  By  considering  dissipation,  however,  log-amplitude, 

phase,  and  log-amplitude-phase  spectra  of  the  scattered  acoustic  wave 
can  be  obtained  over  the  range  of  acoustic  wavelengths  pertinent  to  this 
study.  It  is  shown  that  Tatarskii's  results  for  optical  propagation 
are  good  approximations  for  acoustic  waves.  These  results  are  discussed 
in  Section  A. 

The  equations  governing  acoustic-wave  propagation  are  determined  in 
Section  B.  Assuming  that  variation  of  the  acoustic  refractive-index  is 
much  smaller  than  unity,  which  is  justifiable  under  typical  atmospheric 
conditions,  the  acoustic-wave  equation  and  the  Rytov  technique  developed 
in  Section  B are  applied  in  Chapter  IV  to  describe  the  coherency  of  a 
spherical  acoustic  wave  propagating  vertically  into  the  lower  atmosphere. 


[ A.  Sound  Propagation  in  a Dissipative  Atmosphere 

A first-order  perturbation  solution  of  the  pressure  equation  derived 
I by  Clifford  and  Brown  [1970]  forms  the  basis  for  computing  the  log- 

amplitude,  phase,  and  log-amplitude  phase  spectra  of  the  scattered 
acoustic  wave.  The  atmosphere  is  considered  to  be  dissipative,  and 
tensor  notation  is  used.  The  ground  is  the  reference  level  for  devia- 
tions in  other  physical  parameters  such  as  temperature.  A static 
atmosphere  is  defined  as  one  in  which  ground-level  parameters  exist 
throughout  the  atmosphere. 
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1 . Fluid  Dynamics 

In  the  absence  of  sound,  the  atmospheric  fluid  velocity  is 

”i  ” ^'^i^  + wl , (i  = 1 , 2,  3)  in  the  half  space  x 0;  the  con- 
stant mean-wind  velocity  is  and  w!  is  the  fluctuating  com- 

ponent having  zero  mean.  The  acoustic  pulse  is  assumed  to  be  beamed 
in  the  positive  x-direction  from  the  origin  and,  in  its  presence,  fluid 
motion  is  V^.  = w^.  + u^  where  u^.  is  the  velocity  caused  by  sound 
waves.  Mass  flow  is  expressed  by  pV^. , where  p is  the  fluid  density 
and  Pq  is  the  static  fluid  density.  The  following  notation  is  used: 


9^.  - 9/9x^  , 9^  = 

9/9t, 

and  d^  = d/dt  where  d^  = 

^t" 

V.9..  The 

total  pressure  in 

the  fluid  is  P = p + p ; here 

0 a 

Po 

is  the  constant 

static  pressure  and  p 

a 

is  the  acoustic  pressure. 

In  practice 

'Pa/Pol  I 3nd 

P = 

Pq( 1 + C)  where  j ^1  « 1 . 

The 

shear  strain 

e^j,  fluid  stress 

, and  total  momentum  flux  tensor 

tij  are 

defined  as 

= iO,Vj  t 3.V,) 

(3.1a) 

P. . 

■ (P  - \3kVk)«(j  - Vij 

{3.1b) 

t.  . 
ij 

= "(J  * 

(3.1c) 

In  these  equations,  6..  = 1 for  i = j and  6..  = 0 for  i ^ j, 

1 J 1 J 

and  Ay  and  are  the  volume  and  shear  viscosities,  respectively. 

For  an  ideal  gas  H^,  the  entropy,  is 

= Cy  1 n P - Cp  1 n p ( 3 . 1 d ) 

The  specific  heats  of  air  at  constant  volume  and  pressure  are  Cy  and  Cp 
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respectively,  and  the  Laplace  constant  y is  equal  to  Cp/c^.  The 

• 2 2 
magnitude  of  the  fluid  velocity  c can  be  determined  from  c = y P/P  = 

yR  T , where  R is  the  gas  constant  and  T is  temperature.  If  c 
y ^ y c 0 

and  are  the  reference  fluid- velocity  magnitude  and  temperature, 

2 2 

respectively,  then  c = c^(1+t);  here  t is  the  fractional  variation 
of  temperature  from  the  reference.  In  addition  t = + t' , 

<^t'^  = 0,  and  t'  contributes  to  the  turbulence  fluctuations. 
Assuming  that  the  external  source  terms  are  zero,  the  basic  equations 
of  fluid  motion,  conservation  of  mass,  momentum,  and  adiabatic  propa- 
gation of  sound  are 


9^p  + 9.(pV.)  = 0 

(3.2a) 

9^(PV.)  . a.t..  = 0 

(3.2b) 

dtHe  = 0 

(3.2c) 

Defining  = p/Pq,  = 2v^  + X /p^  and  = P^/CyP  ), 

-ico  t 

then, for  slowly  varying  turbulence,  tt,(x.  ,t)  = e ^ n(x.);  here  ca 

a I 1 a 

is  the  acoustic  angular  frequency.  Combining  the  above  equations  and 
linearizing,  results  in  the  following  equation  for  the  acoustic  pressure 
fiel  d 


2 . 2 

v^n  , . h a , lih  a,a,(„,a,n) 

’ c k^  ’ J'  1 j ' 

a 0 a 


(3.3) 


2 2 

In  this  equation,  V = 9.  is  the  Laplacian  operator,  k, 

1 a 

h = kg  + id,  and  d = kj  n^/(2cQ). 
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. 4.J*  . 


ILj. 


2-  First-Order  Perturbation  Solution  of  the  Pressure  Equation 


Letting  IT  = + 11^  then 


v\  * = 0; 


(3.4a) 


and 


v^n,  . h^n,  . ^ 


(3.4b) 


where  the  fluctuation  in  the  acoustic  refractive  index  n.|  is 


"l  = - 


w 

1 + _>i 

2 Sj 


<"]>  = » 


(3.4c) 


Assuming  this  field  to  be  locally  homogeneous  n^  and  II.|  are  expanded 
in  their  Fourier-Stiel tjes  measure  as 


OO  

i K -y 


— rrr  ’ •y 

= n^(x,0,0)  +yy  e “ “ - 1 


dP(x,K^) 


(3.5a) 


and 


ni{x,y^)  = n^(x,0,0)  + JJ  e - ij  dv(x,K^) 


(3.5b) 


where  n^(x,y^)  real  implies 


dv(x,K^)  = dv  (x,-kJ 


dv(x,-<^)  = dv  (x,K^) 


(3.5c) 

(3.5d) 
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_ dp(x,K  ) 

d(()(x,IC  ) = 77-7 =r-r 

a n (x,y  ) 

o'  '■^a' 


(3.6) 


0 X < 0 

Y(x)  = 1/2  X = 0 


X > 0 


(3.7b) 


a(<  ) = i ( |q|^  + 2k  d)^/^  + i ( hl^  - 2k  d)^/2  (3.7c) 


b(<c^)  = I (lq|^  + 2kgd)''‘'  - (|q|^  - 2k^d) 


1/2  1,,  ,2 


(3.7d) 


Writ!  ng 


Y = Y(<^  - k^) 
+ 'a  a' 


(3.8a) 


Because  |n.,|  « jn^l,  in(n/n^)  n^i/n^.  If  In(n/np)  = x + is^, 

where  x S.|  are  the  perturbations  in  the  log-amplitude  ratio  and 

phase  of  the  scattered  wave,  the  Fourier  Stieltjes  measure  d!ii(x,<  ) of 

a 

the  complex-phase  x + is  related  to  that  of  n^(x,y^)  by 


Let  q = (h^  - then 


q = (a+ib)Y{k^-K2)  + (b+ia)Y(K2-k2) 

Q Ut  Ct  a 


where  Y(-)  is  the  Heaviside  step  function  defined  by 


(3.7a) 


yields 


Y = Y(k^  - <^) 
a a' 


q = (a+ib)Y  + (b+ia)Y^ 


(3.8b) 


(3.8c) 


Using  the  forward-scatter  approximation. 
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A 

cl<J>(x,Kj^)  = ik^y*  dv(x',Kj  e 


(igi+go)(x'-x) 

' ^ dx‘ 


(3.9) 


Here,  = (k^  -bY^  - aYj  and  = (-d  + bY_  + aY^).  Note  that 

2 

a,  b,  g^ , and  g^  are  functions  of 

3 . Log-Amplitude,  Phase,  and  - Ampi i tude-Phase  Jp^c t ra 

Assuming  that  turbulence  is  homogeneous,  the  Fourier-Stiel tjes 
measure  of  the  log-amplitude  '"atio  and  phase  of  the  scattered  acoustic 
wave  da(x,K^)  and  da(x,K^)  are  [Tatarskii,  1961] 


cI4>(x,k^)  = da(x,K^^)  + ido(x,<^) 


(3.10a) 


doi(x,K^)  = 


[d(J)(x,K  ) + d0  (x,-r,J] 


(3.10b) 


do(x,<J  = 


[d<{)(x,r  ) - d(l)  (x,-K  )] 


(3.10c) 


The  log-amplitude,  phase,  and  log-amplitude-phase  spectra  of  the 

acoustic  wave  can  be  obtained  by  computing  certain  auxiliary  spectra. 

From  Chapter  II. A. 3,  the  two-dimensional  refractive-index  spectrum 

F (x' -x" ,<  ) is 
n <r 


<^dv(x',<  )dv  (x",k')^  = 6(k  -k')F  (x'-x",k  )d^K  d^K'  (3.11a) 

\ ' a'  a / A a n a « a 


^d't){x,K  )d())  (x,k')^  = 6(k  -k')F,(x,0,k  )d^K  d^K'  (3.11b) 


2-  J- 

ly  j-i  ty’ 


<(d(t)(x,<Jd<l>(x,K;^))>  = (3.11c) 
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where,  from  Eq.  (3.9) , 


r/  n-^  2 rr  , - [igi(x'-x")  + g„(x'+x"-2x)] 

F, (x.O.K  ) = k:  //  F (x'-x",<  ) e ^ dx'dx" 

I a a JJ  n a' 

(3.12a) 


0 0 


X X 


' -kf  _//■  F„(x'-x",k^)  e 


(igi+go)(x'+x"-2x) 

dx'dx"  (3.12b) 


0 0 


The  above  equations  can  be  evaluated  by  noting  that  F (r,K  ) 

n ^ o 

is  even  in  ^ and,  by  using  the  substitution  x'-x"  = C.  (x'+x")/2  = n. 


.2  X 


" 'if  ® cos(g,0 


and 


-2g2x  / _ g.c 

J dE 

0 


.2  X 


(3.13a) 


_ r "(igi+go)? 

F2(x,0,<J  = (ig^+g^)  f ® 


-2(ig,+g2)x  ^ _ (igi+gj)^ 

J e 1 2 (3.13b) 

0 


Because  the  major  contribution  in  these  expressions  arise  from  ^ _<  L^, 
the  contribution  of  real  terms  in  the  exponent  factor  of  the  integrands 
is  negligible,  and 


(3. 14a) 


L 


f =/  F^(Ji,K:Jcos(g^OdC 


Furthermore,  'jignificant  contribution  to  F^(x,0,k^)  and  F2(x,0,<^) 


results  from  « x,  hence  the  upper  limits  in  the  integrals  in 


Eqs.  (3.13a)  and  (3.13b)  can  be  replaced  by  infinity.  As  shown  in  Chap- 


ter II. A. 3,  the  three-dimensional  refractive-index  spectrum,  ‘l’p('^)»  is 


(JU  » 

2TT<I>n(N/<^+<5)=  / 


(3.14b) 


therefore 


Trk 


F,U,0,<^)  = - 


(3.15a) 


= r#-9?  t-1  * (3.15b) 


Because  F,(x,0,k  ) is  real  and  a function  of  k and  F„(x,0,i.'  ) is 
la  a ^ X 


complex  and  a function  of  k , 

a 


<fd.(x,<  )drx*(x,Kc' )>  = 6(k  -k')F  (x,0,^^)  d^^  d^<' 
<d  (x,K^)da*(x,7;)>  = Mv<;)F5(x,o,7^)  d\  d^K; 

<d.(x,r^)do(x,7-)>  = Mr^+r;)F^5(x,o,Kj  d\  d^^; 


(3.16a) 

(3.16b) 

(3.16c) 


where 


[F^  + Re(F2)J 

^ f 


(3.17a) 
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[Fi  - Re(F2)] 


(3.17b) 


Iin(F2) 


(3.17c) 


Re(F2)  and  Im(F2)  denote  the  real  and  imaginary  part  of  the  complex 

number  F„.  Here  F (x,0,k  ),  F (x,0,k  ),  and  F (x.O.ic  ) are  the 

^ X Oi  s a XS  01 

log-amplitude,  phase,  and  log-amplitude  phase  spectra  of  the  scattered 
acoustic  wave,  respectively,  and  F , F , and  F can  be  calculated 

X s XS 

from  Eqs.  (3.15)  and  (3.17).  Limiting  forms  of  F , F and  F 

X s xs 

are  examined  in  the  following  cases. 


Case  1 : k « k 


In  this  case,  Y = 0,  Y =1,  h k , a = Iql,  k -a 

T “ 3 3 

b~d,  g.  = K^/(2k  ),  and  g ^ 0.  As  a result, 

I cx  a c 


</(2k  ), 
a a 


F,  - 2tt  k^x  <})  (k  ) 
I ana 


(3.18a) 


r -iK^x/k^ 
2 I e « ^ 


^2  = k"  ^n^^) 

^ " (iK^/kJ  " « 

a a 


(3.18b) 


F = 7Tk"x  1 - 
X ^ 


sin(K  x/k,) 
a a 


--O--  - ^Ak) 

(K^X/kJ  " “ 

a a 


(3.19a) 


5 sin(K  x/k,) 

F,  -irklx  1 + , 

" ^ (k  x/k  ) 

a a 


^>n(<  ) 

n'  a 


(3.19b) 


(3.19c) 


\ Q(  / 

The  above  expressions  are  identical  to  those  derived  by  Tatarskii  [1971] 
but  were  obtained  via  a different  approach. 

Case  2:  k = k . 

— a a 

In  this  case,  Y,  + Y =1,  |q|  = (2k, a = b = (k 

g.|  k^,  and  = (k^d)^^^;  therefore. 


F.  = 


j,  -(k  d)^/^x/sinhL(k 
27rkf  e ® I -?-o -l$.(/2  .c  ) 


1 """a 


(k^d) 


1/2  fn' 


(3.20a) 


F2  = 


-(ik^,+(kgd)^/^)2x 

— 4>  (v^  K ) 

ik,  n^  a' 


J 


(3.20b) 


and 


2 -(^d)’/2x 

Fy  = .k^  e ^>^(^<  ) 


sinh((k  d)^'^^x)  -(k,d)^^^x/sin(2k,x)^ 

a _ 0 ” ■ “ 


(kad) 


1/2 


2k. 


(3.21a) 


2 -(k  d)^/2x 
Fs  = Tik^  e 


sinh((k  d)^/^x)  -(k  d)^/^x/sin(2k  x)^ 

a - + e ® ^ 


(kjd) 


1/2 


2k. 


Trk 


/S  ? 


-(k  d)^/^2x 

1 - e cos(k  2x) 


(3.21b) 

(3.21c) 


Case  3: 


< » k . 

a a 


In  this  case  Y=l,  Y =0,  h=k,  a ~ 
a - d = Consequently, 


|q| , = k^,  and 


1 K"  n a 


(3.22a) 


F«  = - ^ ) 

2 < n a 


(3,22b) 


F = 0 
X 


(3.23a) 


F =—$(<) 
s K n'  a' 
a 


(2.32b) 


F = 0 
xs 


(3.23c) 


Normalized  F^,  F^,  and  F^^,  that  is,  = f^^/('^k^x<j>^(<^)) , 

^XS  ^ ^ vs 

1 /2 

normalized  wave  number  k /(k^/x)  in  Fig.  2,  which  shows  the  behav- 

cx  a 

ior  of  these  normalized  spectral  functions.  It  can  be  observed  that 

these  functions  are  oscillatory  at  low  wave  numbers  and  that  for  large 

wave  numbers,  f and  f are  equal  but  f decreases  rapidly.  The 

function  f is  negligibly  affected  by  damping  over  the  entire  range 

of  wave  numbers;  however,  the  shape  of  the  f and  f spectra  is 

X ^ 

altered  slightly. 

These  log-amplitude,  phase,  and  log-amplitude-phase  spectra 
are  largely  dependent  on  wave  numbers  « k^  that  correspond  to 
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Fig.  2.  NORMALIZED  SPECTRAL  FUNCTIONS  VS  NORMALIZED  WAVE  NUMBER, 
k = 10  m-l  and  x = 1000  m, 

d 

turbulent  eddies  which  determine  the  spectrum  of  the  refractive-index 

fluctuations.  Eddies  with  k > k produce  evanescent  waves.  It 

a --  a 

can  be  concluded  that,  because  $ (k  ) k , these  spectra  are 

determined  by  refractive-index  fluctuations  rather  than  by  absorption  in 

the  medium.  As  a result  Tatarskii's  results  derived  for  > « I are 

a 0 

good  approximations  for  L„  > A > I . 

^ 0 a 0 

In  this  section  plane-wave  propagation  was  assumed;  however, 
the  conclusion  that  'the  formulas  derived  on  the  assumption  that 

<<  are  good  approximations  for  are  unaffected,  whether 

the  type  of  propagation  is  plane,  spherical,  or  beam. 


B . Acoustic-Wave  Propagation 

In  this  section,  the  equations  governing  acoustic-wave  propagation 
are  developed. 
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1 

i 

1 


1 . Acoustic  Time- Harmonic  Wave  Equation 

Let  the  acoustic  pressure  impressed  on  the  medium  resulting 
from  an  external  source  be  p (r,t)  = u (r)e  “ . Because  p is 

d d d 

directly  proportional  to  tt,  which,  in  turn,  is  directly  proportional 

a 

to  n,  an  equation  for  u,(7)  can  be  written  similar  to  Eq.  (3.3). 

a 

In  the  lower  atmosphere,  d « k,  and,  as  a result,  the  effect  of 
dissipation  is  neglected  in  this  analysis.  In  practice,  |t|  « 1 and 
iw|/Co  <<  1,  where  |w|  is  the  magnitude  of  atmospheric  fluid  velocity. 
After  preserving  the  terms  linear  in  t and  w^./c^^,  it  can  be  demon- 
strated that  [Tatarskii,  1971] 

V^u  (r)  + k^n^(r)u  (r)=  0 x > 0 (3.24) 

a ad 

where 


n^(r) 

= <(^n^(r)^  + 2n^  (r)6 

(3.25a) 

n^(r)^ 

= ,-<.>  -2  <“'*'> 

\ / c 

(3.25b) 

0 

(r) 

■ --f -<*'iN>) 

(3.25c) 

0 \ / 

<^n^(7))>  = 0 

(3.25d) 

= unit  vector  tangent  to  the  ray 

The  dummy  parameter  6 measures  the  smallness  of  n^(r), 
which  is  the  variation  of  the  acoustic  refractive- index  from  its  average 
value;  generally,  <^n(7)^  is  approximately  unity  and  in.|(r^|«  1. 

The  above  time-harmonic  scalar  wave  equation  describes  the  acoustic- 
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wave  propagation  in  the  half  space  = ((x.y^)  : 0 < x ' 

0 < |y  I < +00} . 

'a' 

When  n(r)  =1,  Eq.  (3.24)  reduces  to  the  Helmholtz  wave  equa- 
tion. Tatarskii  [1971]  demonstrated  that  the  structure  parameters  for 

the  acoustic  refractive  index,  atmospheric  wind,  and  atmospheric  tempera- 
2 2 2 

ture  (C  , C , and  0^)  are  related  as 
n w T 


r2  - ^ 

0 0 


(3.25e) 


This  equation  also  reveals  the  relationship  between  the  structure  param- 
eters of  the  acoustic  refractive  index,  temperature,  and  wind  velocity. 
Eqs.  (3.24)  and  (3.25)  are  applied  in  the  next  section  and  in  Chapter  IV 
to  compute  the  first  and  second  statistical  moments  of  a spherical 
acoustic-wave  propagating  vertically  into  the  atmosphere. 


2.  Rytov's  Technique 


Using  Rytov's  transformation,  u = exp  [Barabanenkov  et  al , 

d 

1971].  The  perturbation  in  complex  phase-amplitude,  4^,  is  computed  as 


Ua(r)  = A^(r)  exp[iS(r)] 


4;(r)  = In  A (r)  + iS(r) 

d 


(3.26a) 


(3.26b) 


When  this  transformation  is  applied  to  the  scalar-wave  equation,  the 


Ricatti  equation  in  4; 


2 2 2 2 

V ip  + (74;)  + k,n  = 0 

d 


(3.27) 
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is  obtained  and  ijj(r)  is  then  expanded  in  power  series  of  a dimension- 
less parameter  6 as 

oo 

i'(r)  = ^ 6^  ip.(V)  (3.28) 

j=0  J 

Using  Eqs.  (3.25a),  (3.27)  and  (3.28)  and  equating  to  zero  terms  of 
like  powers  in  6 result  in 


V^(-o  ''♦o  ■ 

(3.29a) 

- -2k^n, 

(3.29b) 

Note  that  solving  Eq.  (3.29a)  is  equivalent  to  solving 

<'n^'>u  = 0 

ao  a \ / ao 

which  is  Eq.  (3.24)  with  6=0,  and  u^^  = exp  Equation  (3.29b) 

dO  0 

can  be  solved  via  the  technique  outlined  in  Schmeltzer  [1967].  Assum- 
ing that  <(n( r)^  = 1 , 

.Jl(?)  . -2k2/  (3.30) 

V 

ik  |r-r'  I 
a ' ' 

w^iere  ^ (r,r')  = - ^ (3.31) 

47T|?-f| 

'^ao^^^  ” ^ao^*"^  ®xp[''SQ(r) j,  (4>  - i/^q)  measures  the  fluctuation  in 

logarithmic-amplitude  and  phase  of  the  scattered  sound  wave. 
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where  x are  the  perturbations  in  logarithmic  amplitude  and 

phase,  respectively.  Equations  (3.30)  and  (3.31)  will  be  used  to 
determine  the  coherency  of  spherical  acoustic-wave  propagation  in  a 


I 


Chapter  IV 

COHERENCY  DESCRIPTION  OF  ACOUSTIC-WAVE  PROPAGATION 


The  degree  of  interaction  between  electromagnetic  and  acoustic  waves 
depends  on  the  coherency  of  the  latter.  As  a result,  coherency  of  spher- 
ical acoustic-wave  propagation  in  an  inhomogeneous  atmosphere  is  examined 
in  this  chapter  via  a perturbation  theoretic  approach. 

A diagrammatic  technique,  introduced  by  Feynman  [1948],  serves  as 
a basis  for  the  analysis  and  is  used  to  perform  a selective  summation 
of  the  perturbation  series.  Its  basic  philosophy  is  similar  to  that  of 
the  method  of  smoothing  which,  in  turn,  is  related  to  the  Bogoliubov- 
Krylov-Mi tropol ski  method  for  solving  nonlinear  differential  equations 
[Bogoliubov  and  Mitropolski,  1961].  These  techniques  are  reviewed  by 
Frisch  [1968]  and  Nayfeh  [1973]. 

The  coherency  of  an  acoustic  pulse  propagating  in  a turbulent 
medium  is  a measure  of  random  fluctuations  in  and  between  successive 
wavefronts  and  is  determined  by  atmospheric  parameters.  To  simplify 
the  analysis,  it  will  be  assumed  that  the  mean  of  the  acoustic 
refractive  index  <(n[7)^  is  unity  and  that  the  varying  part  of 
this  refractive-index  field  n-jir)  has  a gaussian  probability  density 
function  [Fante,  1975].  The  acoustic  source  is  modeled  as  a delta 
measure  at  the  origin  of  the  coordinate  system.  Source  parameters  are 
restricted  to  be  deterministic. 

A.  First  Moment  of  Acoustic-Wave  Propagation 

The  first  moment  of  a spherical  acoustic  wave  propagating  in  a 
turbulent  atmosphere  is  derived  in  this  section  and  is  then  used  to 
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estimate  the  coherence  length  of  the  wave.  This  length  is  defined  as 

the  altitude  over  which  the  mean  acoustic  field  is  attenuated  by  e'^ 

resulting  from  complex-phase  interference.  Given  the  notation 

- (x,y  ),  the  coherence  length  x is 
a c 

where  G^(x,y  ;0)  is  Green's  function  for  the  spherical  acoustic  wave 

d Ot 

propagating  in  the  turbulent  medium.  To  derive  <(^G  (x,y  ;0)^  , 
Dyson's  equation  [Frisch,  1968J  will  be  solved  by  means  of  Feynman's 
diagrammatic  technique. 


1 . Dyson's  Equation 

Operator  notation  is  used  to  simplify  the  following  presenta- 
tion. Equation  (3.24)  plus  the  source  function  determine  the  amplitude 
c r the  pressure  wave  as 


Here,  S is  assumed  to  be  a delta  measure  at  the  origin  of  the  coor- 

dC 

2 2 2 

dinate  system  and  f h (V  + k n ) is  the  linear  stochastic  differential 

d 

operator.  In  the  above  equation 


f 


(4.3a) 


G 


a 


(4.3b) 


where  P'  and  G'  are  the  varying  components  of  f and  G,.  Using 

d d 

Eqs.  (3.24),  (3.25a),  (4.2),  and  (4.3), 
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I 


(4.4a) 


= 2k; 


(4.4b) 


Therefore,  from  Eqs.  (4.2)  and  (4.3a) 


(<£>  + S’)G,  = S. 


(4.4c) 


G = G S - G f'G 
a ao  ac  ao  a 


(4.4d) 


where  G^^  is  the  inverse  operator  of  It  is  nontrivial  to 

solve  the  above  equation  for  G ; however,  after  formal  iteration, 

d 


G = T'  ;-G ' f 1'^  g'  S 
a j^Q  ^ ao  ^ ao  ac 


(4.5) 


This  is  called  the  Liouville-Neumann  series  and,  because  the  random 
medium  is  semi-infinite,  this  series  diverges  with  finite  probability. 
Because  is  to  be  computed,  a relationship  governing  ^ 

called  the  Dyson  equation,  can  be  obtained  by  the  method  of  smooching; 
therefore,  averaging  Eq.  (4.2)  yields 


(4.6) 


Subtracting  this  result  from  Eq.  (4.2)  yields 


therefore, 


g'  = -g'  (I  .-P  )f'  G'  - g'  f <fG  ^ 
a ao'  1 av'^  a ao  ^ a^ 


(4.7) 


(4.8) 


I 


where  I.  is  the  identity  operator  and  P denotes  the  operation  of 
1 a V 

ensemble  averaging.  By  iteration,  G'  can  be  determined  as 

a 


Substituting  for  G'  in  Eq.  (4.6)  results  in 
a 


<Cg ,)>  = s 

^ a / ac 


(4.10a) 


where 


- 


(4.10b) 


oo 

))1  = X)  I-G'  (I.-P  )s>'p  g' 

\ I ao'  1 av'  f ao  / 


(4.10c) 


Here  HI  is  called  the  mass  operator  and  the  Dyson  operator 
[Frisch,  1968].  Equations  (4.10a)  and  (4.10b)  yield 

The  success  of  Dyson's  method  depends  on  the  fact  that  <Cg  is 
obtained  by  iterating  from  Eq.  (4.11)  after  approximating  ill  instead  of 
iterating  G^  as  in  the  Liouvi 1 le-Neumann  series  and  then  averaging. 
Because  the  expression  for  )ii  is  complicated,  use  of  Feynman's  dia- 
grammatic technique  in  the  following  section  is  justified. 

2 . Series  Solution  of  the  Py^n_  Equation  Using  Feyjiman 's  Diagrams 

Feynman's  diagrammatic  procedure  is  used  to  obtain  a series 
solution  of  the  Dyson  equation.  In  this  technique,  the  source  S is 

dC 

the  Dirac-delta  measure  6(r)  . Noting  that  G^q6(7  - r' ) is  equal  to 
GgQ(r,r'),  the  integral  representation  of  Eq.  (4.2)  then  becomes 


(i„(r,0)  = G^^(r.O)  - 2k^ /g,„( r.r,  )n,  ( r,  )G^( r,  .0)  d’?,  (4.12) 


Formal  iteration  of  this  equation  results  in  the  expansion  of  Eq.  (4.5) 
= '^ao^^’O)  - 2ka/G3o(r,7^)ni  (7^)6^^ (7^,0) 


2 


(4.13) 


The  above  series  can  be  represented  graphically  with  the  following  con- 
ventions : 

C-j  --  Ggg(r.|,r2)  is  denoted  by  a short  line  whose  end 
points  are  r-j  and  r^: 


°ao<''l''‘2>  ' 


1 


^2  - 


2 ' — 

The  random  operator  -f'  = -2k  n,(r)  is  designated 

d I 
y 

by  a dot  placed  on  the  diagram  at  7,  sometimes 
called  a vertex: 


-2k^n,(r) 


--  A dashed  line  joins  the  vertices  for  which 


4k^<ni(r2)ni(ri  )>  is  evaluated: 
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I 


4k^<ni(r2)n^(ri)>  . . . 

•"2  '"l 

C4  - <G,>~ 

For  example, 

ill::; , 

r r^  r'  (4.14) 


The  series  in  Eqs.  (4.5)  and  (4.13)  can  be  expanded  diagrammatical ly  as 


G^(r,0) 


+ 


(4.15a) 


One  of  the  most  significant  attributes  of  this  method  is  its  systematic 
and  explicit  account  of  multiple  scattering.  For  example,  the  third 
term  represents  a wave  excited  at  the  origin;  it  travels  freely  to  and 
scatters  at  the  vertex  r^ , travels  freely  to  and  scatters  at  the 
vertex  r^,  then  travels  freely  to  r.  Averaging  Eg.  (4.15a)  and 
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1 


i 


because  ni(r)  is  a centered  gaussian  random  field,  the  following 
diagrammatic  series  is  obtained  [Tatarskii,  1971]: 

<G,(7,0)>  - + 


+ / 


(4.15b) 


This  series  is  selectively  summed  by  substituting  the  zeroth  term  in 
Eq.  (4.10c)  in  Eq.  (4.11),  and  (4.11)  then  simplifies  to  the  following 
diagrammatic  equation  in  <(^  G 


(4.16) 


On  iteration. 


<Ga(r,0)>  - 


+ . . . (4.17) 

Using  these  equations,  the  coherence  length  of  acoustic-wave  propagation 
is  computed  in  the  next  section,  where  the  condition  of  forward  scatter, 
E^(x),  exists  (because  the  correlation  length  of  the  refractive- 
index  fluctuations  is  on  the  order  of  L^(x)). 

3.  Coherence  Length  of  Propagation 


An  analytic  representation  of  Eq.  (4.16)  is 
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V V 
s s 


<n,{r^)n^(r2)><G^(r^,0)> 


,3~  j3“ 
d r^d 


(4.18) 


Using  the  principle  of  stationary  phase  [Popoulis,  1968]  and  the  mathe- 
matics in  Appendix  B,  this  equation  can  be  solved  for  (r,0)^: 


<G^(r,0)>  = G^^(r,0)e 


— “’2  *\l  '^1  I ^ 


(4.19) 


Here,  the  unperturbed  term  G,„(r,0)  is  damped  by  the  factor 

aO 


sxp[- (x,y^)  I ^/2]  in  a turbulent  medium  because  of  complex- 
phase  interference.  From  Eqs.  (4.1),  (4.19),  and  (A. 20),  the  coherence 
length  is 


54.8  (1.833  - m ) 

d 


C 2 1^2  ^5/3  ^ma 

no  a 0 


6/(1 1-6  m 


(4.20) 


At  altitudes  of  the  acoustic  pulse  x x^,  the  effect  of  turbulence 


is  negligible.  From  Fig.  3,  it  can  be  concluded  that  the  coherence 
length  varies  inversely  with  frequency  and  the  strength  of  turbulence. 
The  incoherence  lengtn  of  propagation  ''j  isdefined  as 


"i 


54.8(1.833  - m ■) 

d 


,6/(ll-6  m^ 


m 


C 2 b5/3  ^ a 

no  a 0 


(4.21) 
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1 


ACOUSTIC  FREQUENCY 
fa  (Hz) 


Fig.  3.  COHERENCE  LENGTH  VS 
ACOUSTIC  FREQUENCY. 

C„o  = Xo  ” ’ ”■ 

Cq  = 340  m/sec,  and  A = 2. 

The  incoherence  length  of  propagation  is  the  altitude  over  which  the 
transverse  coherence  length  of  the  field  is  of  the  order  of  the  inner 
scale  of  turbulence.  For  x^.  <<  x,  the  acoustic  wave  is  incoherent 
as  a result  of  random  fluctuations  in  the  atmosphere.  This  incoherence 
length  of  propagation  is  typically  three  to  five  orders  of  magnitude 
larger  than  the  coherence  length. 

B.  Second  Moment  of  Acoustic-Wave  Propagation 

The  coherence  function  M(r',7")  of  an  acoustic  wave  propagating 
in  a turbulent  atmosphere  is  defined  by 

<(G3(?',0)G*(?',0)>  = ,0)G*q(F",0)M(7' ,?')  (4.22) 
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The  quantity  <^Gg( r' ,0)G^(r" ,0)^  is  the  second  statistical  moment  of 
Green's  function.  This  function  has  been  widely  used  when  a measure 
of  coherence  between  two  points  in  the  same  wavefront  is  required 
[Tatarskii,  1971;  Prokhorov  et  al . , 1975;  Fante,  1975];  however,  the 
definition  and  requirements  in  this  investigation  are  more  general 
because  r'  and  r"  lie  in  the  acoustic  pulse.  To  compute  M(r',r"), 
the  Bethe-Salpeter  equation  [Frisch,  1968]  is  solved  using  Feynman's 
diagrammatic  technique. 


1.  Bethe-Salpeter  Equation 


The  coherence  function  is  a measure  of  fluctuations  between 

two  points  (not  necessarily  in  the  same  plane).  To  compute  M(r',r"), 

the  second  statistical  moment  of  the  acoustic  wave  must  be  determined. 

This  second  moment,  the  ensemble  average  of  the  double  Green's  func- 
★ 

tion  G (r',0)G  (r",0).  This  double  Green's  function  is  defined  as 

d d 


G(r',0)  ® G*{r",0)  = G (r' ,0)G*(r",0) 

d d da 


(4.23) 


that  is,  the  tensor  product  of  Green's  function  and  its  complex  conju- 
gate. In  terms  of  "double  diagrams"  and  Eq.  (4.15a), 


G fr',0)  O G (r",0)  = 

d d 


+ . . . 


(4.24) 
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Each  double  diagram  is  the  tensor  product  of  the  operator  representing 
the  upper  line  and  the  complex  conjugate  of  the  operator  representing 
the  lower  line.  For  example, 


"l  "O  = (.2k2)  //G,,(r,r^)n^(r^)G^^(r^,r-) 


5 Cl  Co 


Vs  Vs 


(4.25) 


The  mean  double  Green's  function  <(^Gg(r',0)  0 G*(r",0)>  is 
expressed  in  the  following  diagrammatic  series: 


<G  (r'  ,0)G*(r",0)> 


where 


(4.26) 


1 ° - (-2kg)  y*y  Gg^(r,r^)GgQ(r^,r^)<(n^(r^)n^(5^)^ 


C Cl  Co 


V V 
s s 


(4.27) 


In  operator  notation,  0 G^  ^ can  be  computed  from  the  following 


Bethe-Sal peter  equation  [Frisch,  1968]: 


<'^i,  ® <>  ' <=a><'^*a>*  <5a><''’*a>'<'^a  ® '^a  > 


where  :i  is  the  "intensity  operator"  because  it  is  used  in  calculating 
the  mean  intensity  of  the  propagating  wave.  Under  a weak  forward- 
scattering condition,  the  diagrammatic  representation  of  this  equation 


<G,  0 G*>  = 


<G^  ® G*>  (4.29) 


The  top  line  denotes  and  the  bottom  one  designates 

This  diagrammatic  equation  can  be  solved  by  iteration  as  follows: 


<G^  0 g;>  = 


(4.30) 


which  is  the  "ladder"  approximation  of  the  solution  to  the  Bethe- 
Salpeter  equation  [Tatarskii,  1971]  . The  Eq.(4.29)  can  be  rewritten  as 

<G^(?',0)G*(7",0)> 

= <G^(?-,0)>  <G*(?’,0)> 


Vs  Vs 


<Ga(ri,0)G*{r2,0)>d^ld^r2 


(4.3i: 


Normalizing  this  equation  by  r ' ,0)G,„(r"  ,0)  and  using  Eq.  (4.22), 

Cl  U o 0 


M(r'  ,r") 


<G^(r'  .0)><G;(r",0)> 

G,„(r'.0)G*  (P',0) 


// 


<G,(P.'-,)><rOp.r,)> 


J J <n,(V,)n,(f,)> 


'■’ao*'-!  .0)G*„(r2.0)H(r,,r3)d=>?,d^2 


which,  used  to  solve  for  M(r',r")  through  Eqs.  (4.19),  (A. 16), 


(4.32) 


(A. 17),  and  (A. 20)  becomes 


M(r',r")  = exp  j-  1 [<(|,j;^(1F')|^)>  + <^|i|^iCr")  | ^) 


' 2^i|;^(r')4;*(r") 


(4.33a) 


Here,  r'  = (x',y^),  r"  = {x"  ,y'^) , and 


I <(ki(r")|^)>  -2<(i|^^(r')il;*(r")^] 

X °° 

4ii^  kf/'/ri-j(Ky„  ~)1'^  ‘1’  (h,<  )dK  dn 

^ J J V oi  a-'am  x ' J an'  a'  a 


(4.33b) 


where 


X = (x'  + x")/2'  and  y = |y'  - y" I 
' -^am  '■'a  -^a' 


2.  Coherence  Function 


The  coherence  function  is  a measure  of  the  sphericity  of  the 
acoustic  wavefronts  propagating  vertically  into  the  turbulent  atmosphere; 
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I 

I 

E- 

1 

I 

I 

t it  also  determines  the  coherency  in  and  between  the  wavefronts  of  an 

I acoustic  pulse.  This  function  is  examined  in  the  following  cases. 

i 

i 

I 

I Case  1 : x « x^ 

r 

[ At  these  altitudes,  the  acoustic  wavefronts  are  almost  coher- 


Case  3:  x^.  <<  x 


At  these  altitudes,  y ■-  -^rfx):  as  a result  k is 

-'am  o''*  a 

integrated  from  0 to  -«r'’^{ri),  and  expanded  in  a 

2 

smal  1 -argument  approximation  as  1 - /4  in  the  right-hand 

side  of  Eg.  (4.33b).  Using  this  expansion  and  Eqs.  (2.23),  (4.21),  and 
(4.33)  the  coherence  function  is 


f 

I M(7',7") 


I 

I 


= exp 


(4.37) 


C . Summary 

Feynman's  diagrammatic  technique  has  been  used  to  describe  mathe- 
matically the  coherency  of  spherical  acoustic  waves  traveling  vertically 
into  the  turbulent  atmosphere.  If  a RASS  is  operating  at  an  acoustic 
frequency  of  1 kHz,  [Brown  and  Keeler,  1975;  Taylor,  1968] 

C ^ = 10  ^ m X = 1 m,  c = 340  m/sec,  A = 2,  and  m^  = 1.33, 

the  coherence  length  can  be  computed  via  Eq.  (4.20)  to  obtain 

g 

x^  = 6.4  X 10  m.  This  result  implies  that  the  coherency  of  an 

acoustic  wave  is  not  affected  by  turbulence  at  x « x when  the 

c 

acoustic  frequencies  are  below  a few  kilocycles  under  typical  atmos- 
pheric conditions. 


1 
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Chapter  V 

SCATTERING  OF  ELECTROMAGNETIC  WAVES  FROM  AN  ACOUSTIC  PULSE 

This  chapter  describes  the  scattering  of  electromagnetic  waves  from 
an  acoustic  pulse  traveling  vertically  into  the  lower  atmosphere. 

Marshall  [1972]  studied  this  phenomenon  in  a static  atmosphere  (without 
turbulence  or  winds).  The  effects  of  turbulence,  mean  wind,  and  tempera- 
ture parameters  are  included  to  obtain  a more  realistic  expression  for 
the  electromagnetic  power  reflected  from  an  acoustic  pulse.  The  funda- 
mentals of  the  physical  phenomenon  are  discussed  in  Section  A,  and  an 
expression  for  the  backscattered  electromagnetic  energy  for  a monostatic 
RASS  is  derived.  In  Sections  B and  C,  the  effects  of  turbulence  and 
mean  atmospheric  parameters  on  the  backscattered  energy  are  considered, 
and  Section  D summarizes  the  results. 

A . c tion  between  Electromagnetic  and  Acoustic  Waves 

A physical  basis  for  the  scattering  of  electromagnetic  energy  from 
an  acoustic  pulse  is  established  in  this  section. 

1 . Maxwell's  Equations 

The  electromagnetic  refractive  index  n^lr)  of  the  medium  is 

defined  by 

ng(r)  = 1 + ng.|(r)  (5.1) 

j 

where  its  varying  part;  typically,  lng^(r)|  <<  1.  The  vari- 

ations in  ng(r)  can  be  caused  by  naturally  occurring  phenomena  in  the 
atmosphere  or  they  may  be  man-made. 
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permeability  is  a constant,  the  electromagnetic  refractive  index  is 


n,(r)  . (5.2) 

where 

e{7)  = Eq  + Ae(r)  (5.3) 

Here,  the  free-space  permittivity  is  a constant,  and  Ae(r)  is  the 

deviation  of  permittivity  e(r)  of  the  medium  from  e . Given  a 

-io)  t ° 

sinusoidal  time- varying  electric  field  E(r)e  ® and  a magnetic  field 
_ _ -iw  t 

H(r)e  , Maxwell's  equations  can  be  expressed  (using  vector  notation) 
as 

V x E(r)  = ipQUgH(r)  (5.4a) 

V x H(r)  = -!£( r)a)gE(r)  (5.4b) 

V • H(r)  = 0 (5.4c) 

V • [E('r)I(r)]  = 0 (5.4d) 

where 

is  the  angular  electromagnetic  frequency  and  f^  is  the  electromagnetic 
frequency.  The  above  equations  can  be  combined  as 

V^E(r)  + (jjgP^e(r)E(r)  + V[E(r)  • V(lne(r))]  = 0 (5.6a) 

The  gradient  term  can  be  neglected  if  is  small  over  one 

wavelength  [Bremmer,  1964].  This  is  a valid  simplification  because 
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this  term,  which  contributes  to  the  depolarization  of  the  electromag- 
netic wave,  is  negligible  under  the  atmospheric  conditions  [Strohbehn, 
1971]  pertinent  to  this  study.  Neglecting  the  gradient  and  applying 
Eq.  (5.3), 


V^E(r)  + a.3yoGoE(r)  = -a)^yQAe(r)E(  r) 


{5.6b) 


2 2 

Using  k = w y c , where  k = 2tt/a  is  the  electromagnetic  wave- 

G 6 0 0 G 6 

number,  the  above  equation  simplifies  to 


V^E(r)  + k^E(r)  = -k^  ^(r) 


(5.7) 


which  can  be  solved  in  a perturbation  series  because  |Ae(T)/GQ|  « 1. 


2.  Perturbation  of  Permittivity  of  Air  by  a Propagating  Sound  Wave 

Neglecting  the  effect  of  water  vapor  [Bean  and  Dutton,  1968], 
the  electromagnetic  refractive  index  above  30  MHz  can  be  approximated  by 


(n  - 1 
e 


(5.8) 


where  pressure  P is  in  Nw/m*^,  temperature  T is  in  °K,  and 

-8  2 

is  a constant  equal  to  77.7  x 10  °K/Nw/m  . Equations  (5.2)  and 


(5.3)  yield 


(c  - 


<"e  - '>  ■ 


(5.9) 


therefore. 
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dc  dP  dT 


e - Gq  P T 


(5.10) 


Because  the  propagation  of  a sound  wave  is  an  adiabatic  process  in  the 
lower  troposphere, 


dP 


dV 


- -Y 


(5.11) 


where  y is  the  ratio  of  specific  heats  nt  air  and  is  its  volume 

density.  Using  the  above  equation  and  the  perfect  gas  equation 
PVd/T  = constant. 


dT 

T 


then 


which,  in  incremental  notation,  is 


ion  and 

the 

/x^j' 

\ dP 

\ ^ > 

1 P 

2K 

re 

dP 

Y 

is 

2K 

re 

AP 

Y 

T 

(5.12) 


(5.13) 


(5.14) 


Defining 


P = p + p 

0 *^3 


(5.15a) 


where  p is  the  static  atmospheric  pressure  and  p is  the  external 

O d 

acoustic  pressure,  then 


AP  = p. 


(5.15bl 
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The  external  source  is  assumed  to  have  a sinusoidal  amplitude  variation 
with  a maximum  value  of  If  P,  is  the  radiated  acoustic 

alTiaX  a 

power  (in  watts),  g,(6)  is  the  gain  function  of  the  acoustic  source, 

d 

and  G is  the  maximum  value  of  g,(0),  then 

dS  d 

ga(0)  = Gas  (5.16a) 

and 

^^(0)  = 1 (5.16b) 

According  to  Landau  and  Lifshitz  [1959],  the  intensity  I(r)  at 
distance  r is 


Hr) 


amax 

2p  c 
0 0 


Airr^ 


(5.17) 


Consequently  combining  Eqs.  (5.15)  and  (5.17)  yield 


(5.18) 


For  the  fol lowing 

K = 77.7  X 10'® 
re 

T = 273  °K,  and 


typical  values  in  the  lower  troposphere, 
°K/Nw/m^,  = 1.23  Kg/m^,  c^  = 330  m/sec, 

Y = 1.4, 


(5.19) 


and 
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(5.20) 


where,  assuming  g (6)  ==  G over  the  effective  scattering  volume  V , 


2 -15 

< = 10  P G 
1 a as 


(5.21) 


3 . Born  Approximation 

Equation  (5.7)  can  be  written  as 


7^E(r)  + k^E(r)  = -k^6g  E(r) 


(5.22) 


for  each  cartesian  component  of  E(r),  where  5^  is  a dimensionless 

dummy  parameter  denoting  a measure  of  deviation  of  relative  permittivity 

from  unity.  The  electric-field  component  E(r)  can  be  expanded  in 

a power  series  of  6 as 

e 


E(r)  = E <5^  E.(r) 

j=o  ^ 


(5.23) 


which,  when  substituting  into  Eq.  (5.22)  and  equating  terms  of  the  same 


order  in  6^,  results  in 


2 2 
V E„  + kfE„  = 0 
0 e 0 


(5.24) 


V^E,  + k^E,  + k^  c = 0 

1 el  e \ Eq  / 0 


(5.25) 


m e m e \ e I m-1 


(5.26) 


In  the  above  set  of  equations,  E-j  produces  only  single-scattering  and 
higher  order  terms  yield  multiple-scattering  effects;  m^^  order 
scattering  is  a result  of  E.,  m ^ j ^ 0.  A sufficient  condition  for 
the  validity  of  the  above  series  is 


1 — d^  <1 

2 Gq  e scv 

max 


where  d^^^  is  the  diameter  of  the  scattering  volume.  If 


1 tel  « 1 

2 e I e scv 

max 


then , 


E(r)  . E^(r)  + E.,(r) 


(5.27) 


This  is  known  as  the  Born  approximation  [Bremmer,  1964].  If  is 
the  transmitted  electromagnetic  power  (in  watts),  is  the  maximum 

gain  of  the  electromagnetic  transmitter,  and  rig  is  the  free-space 
impedance,  then 


where 


_ i(k„r-u  t) 

EJO  • E„„(f)  e ' ® 


4<e)  . E* 


(5.28) 


2%  47Tr2 


(5.29) 


Using  Green's  function  theory  and  imposing  the  radiation  condition  [Yeh 
and  Liu,  1972]  the  scattered  field  at  the  receiver  for  monostatic 
geometry  of  the  RASS  is 


I 


traveling  in  a turbulent  atmosphere  is  computed  for  the  RASS  mono- 
static geometry  illustrated  in  Fig.  5.  The  half-power  beamwidths  of 
the  acoustic  source  and  the  electromagnetic  antenna  are  9 and  9 , 

3 0 


Fig.  5.  RASS  MONOSTATIC  GEOMETRY. 
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i 

\ 


i 


I 


respectively,  and  0 is  the  minimum  of  each.  The  effective  scattering 
volume  V is  limited  by  two  spheres  of  radii  R,  and  R_  centered 
at  the  origin  and  by  a cone  of  vertical  angle  6.  If  n is  the 

d C 

number  of  cycles  in  the  acoustic  pulse,  then 


Rn  - Ro  = n 
1 t ac  a 


(5.31) 


The  maximum  gains  of  the  electromagnetic  transmitting  and 
receiving  antennas  are  and  G^,  respectively.  These  antenna 
patterns  are 


g^(9)  - 

fj(e) 

(5.32a) 

^^(0)  = 1 

(5.32b) 

f?(e) 

(5.32c) 

- 1 

(5.32d) 

Here  g..,(6),  g,.(0)>  and  the  acoustic  source  pattern  g,(e)  defined  in 

I L a 

Eqs.  (5.16)  are  considered  to  be  substantially  constant  over  V , as 
is  E^^(r)  from  Eq.  (5.29).  These  simplifying  assumptions  do  not 
significantly  affect  the  calculations  because  the  variation  in  these 
patterns  and  amplitude  over  is  small.  If  the  effective  area  of 

the  receiving  antenna  is  and  the  received  power  is  P^, 


and 


A = — G 


(5.33a) 


i 


1 


I 
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A 2ri„ 

r v 


(5.33b) 


The  diameter  of  the  scattering  volume  is 


d = max  i R0  , n A z 
scv  I ac  a ) 


(5.34) 


where  R x is  the  height  of  the  acoustic  pulse.  As  a result,  it  can 
be  noted  from  the  previous  subsection  that  the  condition  for  the  validity 
of  the  Born  approximation  is 


K^kgX0^/2  « 1 


if  d = x0 
scv 


“scv  • "ac^a 


These  conditions  are  typically  valid  in  the  lower  troposphere. 


B . Received  Power  in  the  Presence  of  Turbulence 

In  computing  the  scattered  electromagnetic  energy  from  an  acoustic 
pulse  in  a turbulent  atmosphere,  the  mean  atmospheric  parameters  are 
assumed  to  be  absent.  The  effect  of  these  mean  parameters,  however,  is 
considered  in  the  next  section.  From  the  results  obtained  in  Eqs. 
(3.32)  and  (5.20) 


= — Re  / e 

r I 


Lik  r+,jj(r)-m  (r)]' 


(5.35) 


Based  on  this  equation  plus  (4.22),  (4.33a),  (5.28),  and  (5.30)  and 
imposing  the  Bragg  condition  2kg  = k^. 
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valid  even  in  the  presence  of  atmospheric  turbulence. 


Case  2:  x <<  x <<  x. 

c 1 

In  this  range,  6/2,  can  be  approximated  by  ird^  in 

the  integrand  of  Eq.  (C.9);  then, 

2d 

I(x,  6/2,  y ) = (2Tr)(TTd^)  / y M(x,y  )dy  (5.42) 

' ■'am  a’ J ■'am  ' ■'am'  ■'am 

0 

Because,  generally,  2d  »'^-(x),  the  upper  limit  of  this  integral 

d 0 

can  be  extended  to  infinity.  Based  on  Eq.  (4.36), 


r 


\ — 


(8/3  rn  ^ ^ 

iT/6  - 


x\(l1/6-m  1(6/5) 


and 


I(x,  e/2,  = (Trd2)2 


am'  ' "a' 

max 


As  a result. 


I(x,  e/2,  y_) 


am' 


Kx,  6/2.  yj 


= 0.226 


max 


/ 8'^  - "af 
^Tw^;  \Tr/ 


(^) 


(ll-6m^)/5 


(5.43) 


(5.44) 


(5.45) 


^rmax  (5.41)  is  the  maximum  received  power  in  the 

absence  of  turbulence,  then  from  Eqs.  (5.39)  and  (5.45),  it  follows  that 


W 

P 

rmax  I(x,  0/2,  y^^^) 


max 


0.226 


(ii/6  - »,j  y d,  j 


/x^|’'-6"a)/5 


(5.46) 


from  which  it  can  be  observed  (Fig.  6)  that  the  received  power  decreases 
as  altitude  increases. 
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Fig.  6.  NORWLIZED  RECEIVED  POWER  VS  ALTITUDE  OF  ACOUSTIC 


PULSE.  , C2  = 10-6  m-2/3. 


f,  = 10^  Hz,  A = 2, 


1 m, 


' '0  ....  5 

= 0.1  rad,  X « x <<  x . 

c 1 


340  m/sec. 


Case  3:  x^  <<  x 

At  these  altitudes,  it  can  be  shown  by  using  Eqs.  (4.37),  (5.39), 
(C.8)  and  (C.9)  that 


rmax 


(ll/6-inJ 

d 


:5.47) 


It  can  be  concluded  from  the  results  obtained  for  the  three  cases 
considered  above  that,  in  the  lower  troposphere  (x  ~ 10  km),  the 
effect  of  turbulence  on  the  received  RF  power  in  a RASS  is  negligible 
at  acoustic  frequencies  less  than  a few  kilohertz.  At  frequencies 
greater  than  i 10  kHz,  however,  turbulence  must  be  considered  when 
evaluating  the  RASS  performance  when  the  altitude  of  the  acoustic  pulse 
is  greater  than  about  1 km. 
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C . Interaction  between  Electromagnetic  and  Acoustic  Waves  in  the 
Presence  of  Mean  Atmospheric  Parameters 


In  addition  to  the  fluctuating  parameters,  the  atmosphere  is 
characterized  by  such  mean  parameters  as  temperature  and  winds.  In 
the  lower  atmosphere,  a linear  temperature  profile  can  be  observed  at 
heights  up  to  the  first  10  km;  its  temperature  gradient,  typically 
-6.5  °K/km  [Yeh  and  Liu,  1972],  tends  to  disperse  the  acoustic  wave 
train.  The  atmosphere  is  also  characterized  by  vertical  winds  that 
alter  the  effective  velocity  of  sound  and  by  horizontal  winds  that  shift 
the  focus  of  electromagnetic  energy  from  its  origin.  In  the  analysis  to 
follow,  these  winds  are  assumed  to  be  constant,  and  each  parameter  is 
studied  separately. 

1 . Received  Power  in  an  Atmosphere  with  a Linear  Temperature 
Profi le  ~ ~ 

The  linear  temperature  profile  in  the  lower  atmosphere  can 
be  described  as 


T(r)  = T'  + a^r  (5.48) 

where  a^  is  the  temperature  gradient,  T'  is  the  temperature  at 
ground,  and  T(r)  is  the  temperature  at  r.  The  beamwidth  of  the 
acoustic  source  is  assumed  to  be  sufficiently  narrow  so  that  the 
temperature  variations  across  the  acoustic  wavefront  can  be  neglected. 
This  is  a valid  assumption  for  acoustic-wave  propagation  in  the  tropo- 
sphere [Yeh  and  Liu,  1972]  and  for  beanwidths  less  than  1 radian.  The 
endpoint  of  the  vector  r falls  within  the  half-power  beamwidth  of  the 
acoustic  source.  Because  I a r|  <<  T'  at  altitudes  below  a few 

I g I 


kilometers,  T'  will  be  replaced  by  the  average  temperature  in  the 


ensuing  calculations. 


From  Tatarskii  [1971]  and  Eqs.  (3.24),  (3.25)  and  (5.48),  the 


phase  S.p(r)  of  the  acoustic  wave  is  given  by 


ST(r)  = k^r  - r\r 


(5.49a) 


where 


a k^ 

= 


(5.49b) 


In  the  geometry  in  Fig.  5 (page  58),  the  strength  of  the  RASS  received 


signal  is  proportional  to  the  magnitude  of  F,  where 


I 

F = y"  cos  S.|-(r)  e 


-i(2k^r-a) 


(5.50) 


and  a is  a constant  introduced  to  simplify  the  calculations  of  |F|  . 


In  the  absence  of  temperature  gradients. 


(5.51) 


For  a linear  temperature  profile  in  the  lower  atmosphere,  F calculated 


in  Appendix  D is 


t s(n(nR„+a)tSf(AJ  - 5f(BJ] 


66 


+ cos(nRp-a)[C^(Ap)  - C^(Bp)] 

+ sin(nR^-a)[S^(Ap)  - S^(Bp)] 

+ i cos(riR^+a)[-S^(A^)  + S^(B^)] 

+ sin(nR^+a)[C^(A^)  - C^(B^)] 

+ cos(nRp-a)[S^(Ap)  - S^(Bp)] 

- sin(nRp-a)[C^(Ap)  - C^(Bp)]  j (5.52) 

where  C^(-).  S^(-)  are  Fresnel  cosine  and  sine  integrals,  the  nor- 
malized temperature  gradient  B is  defined  as 


and 


B = 


(5.53) 

(5.54a) 

(5.54b) 

(5.55a) 

(5.55b) 

(5.55c) 

(5.55d) 


(5.56a) 


(5.56b) 


(5.56c) 


(5.56d) 


? 2 

Defining  the  power-reflection  coefficient  to  be  IFl  /(n,^X,/2)  , 
it  can  be  observed  that  a match  between  the  electromagnetic  and  acoustic 
wavelengths  depends  on  the  temperature  difference  over  the  length  of  the 
acoustic  wave  train.  For  an  acoustic  pulse  at  an  altitude  x,  the 
reflection  of  electromagnetic  energy  is  maximized  when 


This  represents  a modified  Pragg-s^atter  condition  when  a RASS  is  operat- 
ing in  the  presence  of  a linear  temperature  profile.  For  |6|  « 1, 
under  modified  Bragg-scatter  condition  and  small  wave-train  lengths,  the 
power-reflection  coefficient  is  almost  unity.  It  can  be  seen  in  Fig.  7 
that  the  wavelength  ratios  that  produce  a peak  in  the  power-reflection 
coefficient  for  positive  and  negative  temperature  gradients  are  differ- 
ent; this  can  also  be  observed  in  Eq.  (5.57).  The  conclusion  drawn  from 
these  results  is  that,  ^or  short  wave-train  lengths  and  wavelengths 
below  a few  meters,  typical  temperature  gradients  in  the  lower  atmos- 
phere have  a negligible  effect  on  the  received  power  of  a RASS. 


> 
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LOSS  (dB)  LOSS  (dB) 


a.  Negative  temperature  gradient 


X,/Xg  RATIO  OF  electromagnetic  TO  ACOUSTIC  WAVELENGTH 


b.  Positive  temperature  gradient 


Fig.  7.  EFFECT  OF  LINEAR  TEMPERATURE  PROFILE  ON 
POWER-REFLECTION  COEFFICIENT  VS  WAVELENGTH  RATIO; 
altitude  = 1000  X,,  pulse  length  = 40  X . 

a d 


2 . Received  Power  in  an  Atmosphere  with  a Linear  Temperature 
Profile  and  Mean  Vertical  Wind 

In  the  presence  of  a linear  temperature  profile  described  by 
Eq.  (5.48)  and  a constant  mean  vertical  wind  the  phase  of  an 

acoustic  wave  propagating  in  the  lower  troposphere  can  be  determined 
from  Eqs.  (3.24),  (3.25)  and  (5.49b).  This  phase  is 


S,(r)  = k^r 


('  - ¥)  - 


nr 


(5.58) 


In  the  geometry  in  Fig.  5 (page  58),  the  strength  of  the  RASS  received 
signal  is  proportional  to  the  magnitude  of  F^,  where 


-f 


cos  S (r)  e 

V 


-i(2k  r-a) 


dr 


(5.59) 


and  a is  a constant.  In  the  absence  of  temperature  gradients  and  winds. 


Defining 


<“«> 


a = 

"o 


(5.60) 


(5.61) 


then,  F^  can  be  calculated  as  in  Appendix  D, 


* s1n(nR^/  »):Sf(A„,)  - 
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Figure  8 illustrates  the  power-reflection  coefficient  defined  by 
2 2 

^ function  of  whose  shape  is  determined 

by  the  temperature  gradient.  The  optimal  ratio  of  A /A  , however, 

6 8 

is  obtained  by  the  temperature  gradient,  altitude  of  the  acoustic 
pulse,  mean  temperature  T , and  normalized  mean  vertical  wind  a . 

0 V 

The  modified  Bragg-scatter  condition  is 
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XgAfl  O'"  electromagnetic  to  acoustic  wavelength 


Fig.  8.  EFFECT  OF  LINEAR  TEMPERATURE  PROFILE  AND 
VERTICAL  WIND  ON  POWER-REFLECTION  COEFFICIENT  VS 
WAVELENGTH  RATIO:  Altitude  = 1000  A , pulse 
length  = 40  A . ^ 
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When  ^e  ” ^ condition  for 

the  maximum  return  of  electromagnetic  energy  in  a static  atmosphere. 

It  can  be  concluded,  as  in  the  previous  subsection,  that  the  loss  in 
received  electromagnetic  power  is  negligible  in  the  presence  of  tempera- 
ture gradients  and  mean  vertical  winds. 

3 . Received  Power  in  the  Presence  of  Horizontal  Winds 

Horizontal  winds  displace  the  acoustic  wavefront  laterally. 

/ 

This  results  in  a shifting  of  the  focus  of  reflected  electromagnetic 
energy  from  its  origin,  thereby  resulting  in  a decrease  in  the  received 
signal  level  when  a monostatic  RASS  is  used.  The  effect  of  constant 
mean  horizontal  winds  w^  is  investigated  as  follows. 

Consider  an  electromagnetic  wavefront  interacting  with  an 
acoustic  wavefront.  The  phase-matching  error  between  the  two  wavefronts 
is  illustrated  in  Fig.  9.  The  error  between  paths  P and  P'  with 
angular  separation  0 is  2(R-R'),  where  R and  R'  are  the  vectors 
corresponding  to  the  two  paths.  The  factor  2 is  the  result  of  the 
forward  and  reverse  travel  of  electromagnetic  energy.  Because  the  phase- 
match  error  is  small  for  small  9,  the  electromagnetic  energies  from 
P and  P'  add  nearly  in  phase  at  the  receiver  located  at  the  origin 
of  the  coordinate  system.  If  P is  fixed  and  0 is  increased,  the 
phase-match  error  will  become  larger,  and  this  will  increase  the  phase 
difference  in  the  electromagnetic  energies  reflected  from  the  two  paths. 
When  this  error  is  half  of  an  electrical  wavelength,  the  energies  will 


I 


T.R.S 
K— H 

DISPLACEMENT 


Fig.  9.  SHIFT  OF  ACOUSTIC  WAVEFRONTS  IN  THE  PRESENCE  OF 
HORIZONTAL  WINDS.  T - Radar  Transmitter,  R - Radar 
Receiver,  S - Acoustic  Source. 

cancel.  If  9 is  further  increased  this  error  may  become  a full 
electrical  wavelength  and  the  energies  will  add  constructively. 

From  this  discussion,  it  can  be  concluded  that  a smaller  beamwidth  b 
reduces  the  susceptibility  of  reflected  electromagnetic  energy  to 
horizontal  winds. 

Each  acoustic  wavefront  in  an  acoustic  pulse  contributes  a 

phase-match  error.  This  error  is  substantially  constant  over  an 

acoustic  pulse  containing  a few  tens  of  cycles  for  typical  values  of 

beamwidths  and  horizontal  wind  velocities;  however,  the  alignment 

between  electromagnetic  and  acoustic  wavefronts  is  important. 

From  Tatarskii  [1971],  and  Eqs.  (3.24)  and  (3.25),  the  phase 

S (?)  of  the  acoustic  wave  is 
w ' 
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(5.66) 


S„(r)  = kjr(l 


where  r = T/r  and  \w/  is  the  mean  wind  vector.  For  the  geometry 
in  Fig.  5 (page  58),  the  magnitude  of  the  RASS  received  signal  is  pro- 
portional to  F^,  where 


/ 


cos 


S (r-) 
w ' 


-i(2k  r'-a)  3_ 
e ® d>' 


(5.67) 


as 


and  a is  a constant.  In  the  absence  of  horizontal  winds  w,  , the 

h 

maximum  value  of  F is 
w 


F 


wmax 


2 

TTX  n 


A 

ac  a 


1 

2 


(5.68) 


where  x is  the  height  of  the  acoustic  pulse.  Note  that,  for  x ->>  A 

d 

and  w^/c^  <<  1,  the  acoustic  wavefronts  shift  laterally  through 
kaXWh/Cp  radians  in  the  presence  of  horizontal  winds.  The  effects  of 
horizontal  winds  on  the  received  power  are  negligible  when 

Ae  = 2Ag,  then  from  Eq.(E.15)  in  Appendix  E,  for  small  pulse  length, 

P 20.37  cos^[k  x0w./(2c  ) - 3tt/4] 

— = f — (5.69 

rmax 


Here  P^,  and  P are  the  power  received  in  the  presence  of  hori- 

rw  rmax 

zontal  winds  and  in  a static  atmosphere,  respectively.  The  functional 
behavior  of  Eq.  (5.69)  is  illustrated  in  Fig.  10. 
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Fig.  10.  ASYMPTOTIC  NORMALIZED  RECEIVED  POWER  VS 
NORMALIZED  ALTITUDE  FOR  ACOUSTIC  PULSE  ALTITUDES 
LESS  THAN  COHERENCE  LENGTH  IN  THE  PRESENCE  OF 
CONSTANT  HORIZONTAL  WINDS.  0 = 0.5  rad. 


If  X = 200  X,,  w,/c„  = 0.01,  and  0 = 0.5  rad,  then  P /P 

a h 0 rw  rniax 

4 X 10”^.  When  X,  = 4 m and  c„  = 340  m/s,  then  f = 85  Hz  which 

d 0 0 

is  the  operating  frequency  of  the  Stanford  RASS  [Frankel  and  Peterson, 


1976];  the  loss  in  received  power  as  a result  of  horizontal  winds  at 

_9 

an  altitude  of  800  m is  P /P  ~ 4 x 10  . If  the  operating 

rw  rmax  ^ 
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and  if  the  beamwidth  remains  constant  at  6 = 0.5  rad,  then  at 

X = 200  ^ = 68  m,  P /P  = 4 x 10‘^. 

a rw  rmax 

The  above  decrease  in  the  received  signal  level  results  because 
horizontal  winds  shift  the  focal  point  of  the  reflected  electromagnetic 
energy  from  the  origin,  and  this  shift  is  a measure  of  the  horizontal 
wind.  Applying  the  principle  of  specular  reflection  (Fig.  11),  a 
bistatic  radar  configuration  with  space-diversified  receiving  antennas 
can  be  used  to  collect  this  reflected  energy  [Frankel  et  al . , 1977] 
and  measure  horizontal  winds. 


Fig.  11.  SPECULAR  REFLECTION  OF  ELECTROMAGNETIC  ENERGY  FROM 
ACOUSTIC  WAVEFRONT.  T - radar  transmitter,  R - radar  receiver 
at  the  origin,  R^  - displaced  receiver,  S - acoustic  source, 
x - altitude  of  the  acoustic  wavefront. 
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D . Summary 


Scattering  of  electromagnetic  waves  from  an  acoustic  pulse 
traveling  in  a realistic  atmosphere  has  been  considered.  The  effects 
of  such  atmospheric  parameters  as  turbulence  and  mean  winds  were  con- 
sidered independently.  These  phenomena  can  be  combined  to  analyze 
their  effect  on  the  scattering  of  electromagnetic  waves;  however,  the 
basic  physical  ideas  remain  unchanged.  The  asymptotic  results  of  the 
different  cases  can  be  combined  to  obtain  a composite  view,  but  only  at 
the  expense  of  clarity. 

It  has  been  concluded  that  turbulence  has  little  influence  on  the 
scattered  electromagnetic  energy  of  acoustic  pulses  with  carrier  fre- 
quencies below  a few  kilocycles  and  at  heights  up  to  a few  kilometers. 
Mean  vertical  winds  and  temperature  gradients  are  functions  of  altitude, 
and  their  effects  on  the  RASS  signal  can  be  minimized  by  establishing 
a modified  Bragg-scatter  condition. 

Horizontal  winds  displace  the  acoustic  pulse  downwind  resulting  in 
received  signal  levels  below  those  computed  in  Eq.  (5.41).  By  using  a 
bistatic  RASS  configuration,  however,  the  principle  of  specular  reflec- 
tion of  RF  energy  from  the  acoustic  wavefronts  can  be  applied  to  measure 
horizontal  winds.  In  this  RASS  geometry,  the  radar  antennas  and  acoustic 
sources  are  aligned  in  the  wind  direction  --  the  transmit  antenna  unwind 
from  the  acoustic  source  and  the  receive  antenna  downwind.  If  the  sepa- 
ration between  the  antennas  is  much  less  than  the  altitude  of  the  acous- 
tic pulse,  the  received  signal  level  can  be  obtained  from  Eq.  (5.41). 
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Chapter  VI 

DOPPLER  SHIFT  IN  THE  REFLECTED  ELECTROMAGNETIC  SIGNAL 

The  electromagnetic  energy  reflected  from  a traveling  acoustic  wave 
train  shifts  from  the  frequency  of  the  incident  electromagnetic  wave. 

This  is  called  the  doppler  shift,  and  it  depends  on  atmospheric  param- 
eters and  on  the  position  of  the  acoustic  wave  train  with  respect  to 
the  RASS  geometry. 

The  influence  of  mean  winds  on  the  received  doppler  signal  will  be 
discussed  in  Section  A.  In  Section  B,  the  doppler  shift  for  a bistatic 
RASS  is  described.  Because  deviation  from  unity  in  the  acoustic  refrac- 
tive index  caused  by  turbulence  is  smaller  than  that  resulting  from 
typical  mean  wind  fields  in  the  troposphere  (by  at  least  one  order  of 
magnitude),  the  effect  of  turbulence  is  ignored  in  this  chapter. 

A . RASS  Doppler  Measurements  in  the  Presence  of  Mean  Wind  Field 

In  monostatic  geometry,  the  doppler  shift  in  the  received  electro- 
magnetic signal  Af(x)  from  an  acoustic  signal  at  an  altitude  x in 
the  absence  of  winds  is 

1 

I 

(6.1) 

1 
! 

c(x)  = K^[T(x)]^/2  (6.2) 

where  = 20.053  [North,  1974J  and  the  speed  of  sound  c(x)  at  x 
is  measured  in  meters  per  second.  From  Eqs.  (6.1)  and  (6.2),  the 
temperature  T(x)  at  x can  be  measured  by  determining  Af(x). 


A .r  / . . \ 

ai  V a;  - 


'-V  A/ 
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Because  horizontal  winds  contribute  negligible  error  in  the  compu- 
tation of  T(x)  for  a monostatic  RASS  geometry,  only  vertical  winds 
w^  need  be  considered.  The  doppler  frequency  in  the  presence  of  these 
winds  is 


Af(x)  = - [c(x)  + w ] (6.3) 

^e  ^ 

Knowledge  of  w^  and  Af(x)  produces  a correct  indication  of  tempera- 
ture. If  vertical  winds  are  ignored,  the  error  in  the  calculation  of 
T(x)  is  [North,  1974]  approximately 

AT(x)  = 1-7  (6.4) 

and  this  error  becomes  significant  when  the  vertical  wind  is  greater 
than  1/2  m/sec. 

The  functional  error  in  the  computation  of  the  height  of  the  acous- 
tic pulse  will  be  on  the  order  of  w^/c^.  This  error  should  also  be 
corrected  because  it  becomes  increasingly  significant  as  the  magnitude 
of  the  vertical  wind  and  the  altitudes  grow  larger.  It  can  be  elimi- 
nated, however,  by  a suitable  averaging  process,  one  of  which  is 
described  by  North  [1974]. 

In  the  bistatic  geometry  of  Fig.  12,  the  doppler  shift  Af(x)  in 
the  received  electromagnetic  signal  caused  by  the  traveling  acoustic 
signal  is 

Af(x)  = - c(x)  cos  [u(x)]  (6.5) 
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Fig.  12.  DOPPLER-TRACKING  BISTATIC 
RASS  GEOMETRY.  T - transmitter, 

S - acoustic  source,  R - receiver, 
X - altitude  of  the  acoustic  wave- 
front,  d,  = TS  = SR. 

d 


which  differs  from  Eq.  (6.1)  because  of  the  additional  multiplicative 
factor  cos  [a(x)],  where 


cos  [a(x)] 


(6.6) 


and  d is  the  separation  between  the  acoustic  source  and  the  electro- 

d 

magnetic  transmitter  or  receiver.  At  an  altitude  x,  the  temperature 
T(x)  can  be  computed  via  Eqs.  (6.2),  (6.5)  and  (6.6).  At  x d , 

d 

cos  [o(x)]  =;  1 and  temperature  can  be  determined  from  Eqs.  (6.1)  ar-? 

(6.2). 

In  the  presence  of  horizontal  winds,  w^,  the  doppler  ► • 
the  bistatic  RASS  for  small  values  of  a is  given  by 
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Af(x)  = - jp  c'(x)  cos  [a{x)J 
''e 

where 

c'(x)  = c(x)  cos  sin  (6.7b) 

In  Eq.  (6.7b)  = ^h^^o  typically  much  smaller  than  unity. 

Hence  c'(x)  = c(x),  and  therefore  it  can  be  concluded  that  the  effect 
of  horizontal  winds  on  the  doppler  shift  in  a bistatic  RASS  is  insignif- 
icant. 

B.  Averaged  Doppler  Shift  of  the  RASS  in  a Bistatic  Con)  wjuration 

A static  atmosphere  is  assumed  in  this  analysis.  Tlie  histatic 
geometry  of  Fig.  12  is  considered.  As  the  acoustic  pulse  spreads  over 
a distance  of  n^^x,  in  the  vertical  direction,  the  effective  altitude 
Xp  of  the  acoustic  pulse  should  correspond  to  the  maximum  intensity  in 
the  doppler-frequency  spectrum  [Gill,  1965].  Note  that  ’‘2  - *p  ' *1’ 
where  x^  and  are  the  heights  of  the  upper  and  lower  edges  of  the 
acoustic-pulse  train.  Beam  power  patterns  of  the  electromagnetic  trans- 
mitter and  receiver  antennas  detennine  x^.  For  example,  if  this 
pattern  has  a gaussian  variation,  Xp  will  be  close  to  the  top  of  its 
edge  x^  at  low  altitudes  of  the  acoustic  pulse.  At  higher  altitudes, 
however,  variation  in  the  beam  power  pattern  over  the  spread  of  the 
acoustic  wave  train  in  the  vertical  direction  is  negligible;  therefore, 

Xp  = (x^  + x^)/2. 

If  the  doppler  shift  is  obtained  by  an  averaging  measurement  of 
■f(x)  when  the  acoustic  pulse  occupies  two  different  positions  in  the 
vertical  direction,  then 

i 
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where  heights  x^i  and  indicate  the  maximum  power  return  of 

electromagnetic  energy  from  the  acoustic  pulse  at  the  two  positions. 
Note  that  X is  any  point  between  x^^  and  x^^  and  that  the  speed 
of  sound  c(x)  is  assumed  to  be  constant  over  Xp^  to  Xp^.  At  high 
altitudes,  cos  [a(x)]  = 1 and  the  averaged  doppler  shift  can  once 
again  be  computed  by  Eq.  (6.1).  At  low  altitudes  and  for  a finite 
spread  of  the  acoustic-wave  train,  however,  Eq.  (6.8)  should  be  used. 
Based  on  Eqs.  (6.6)  and  (6.8),  Af,.,(x)  can  be  simplified  to 

aV 


Af  (x) 
av  ' 


2c(x) 

VvrV 


(6.9) 


Consequently,  the  temperature  at  x can  be  computed  from  Eqs.  (6.2)  and 
(6.9). 
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Chapter  VII 
CONCLUSION 

In  this  stu(iy,  the  interaction  of  electromagnetic  and  acoustic 
waves  in  a stochastic  atmosphere  has  been  discussed  for  a monostatic 
RASS  geometry.  Spectral  expansions  for  the  physical  parameters 
describing  acoustic-wave  propagation  in  a dissipative  atmosphere  were 
developed  from  basic  conservation  equations.  Using  Rytov's  transforma- 
tion and  Feynman's  diagrammatic  method,  it  has  been  theoretically  demon- 
strated that,  for  acoustic  pulses  with  carrier  frequencies  below 
a few  kilohertz  and  propagating  under  typical  atmospheric  conditions, 
turbulence  has  little  effect  on  the  strength  of  the  backscattered  RASS 
radio  signal  at  heights  up  to  a few  kilometers.  This  result  implies 

that  focusing  of  radio- frequency  energy  by  the  acoustic  wavefronts  is 

-2 

primarily  a function  of  sound  intensity  which  decreases  as  x , where 
X is  the  altitude. 

The  effect  of  mean  vertical  winds  and  temperature  profiles  on  the 
strength  of  the  received  signal  was  also  shown  to  be  insignificant. 
Horizontal  winds  shift  the  focus  of  the  reflected  electromagnetic  energy 
from  its  origin,  thereby  resulting  in  a decrease  in  the  received  signal 
level  when  a monostatic  RASS  is  used.  The  principle  of  specular  reflec- 
tions, however,  can  be  utilized  in  a bistatic  radar  configuration  with 
space-diversified  receiving  antennas  to  measure  remotely  the  horizontal 

wind  component  in  the  lower  troposphere. 

The  doppler  performance  of  the  RASS  has  also  been  described.  It  has 

been  observed  that  the  doppler  shift  is  sensitive  to  the  presence  of 
winds,  but  this  sensitivity  can  be  eliminated  by  an  averaging  process. 


» 


A modified  expression  for  the  doppler  shift  was  also  obtained  for  a 
bistatic  RASS  geometry,  which  takes  into  account  the  finite  spread  of 
the  acoustic  pulse  in  the  vertical  direction. 

This  research  has  been  theoretical  in  nature.  The  conclusions 
drawn,  however,  need  to  be  verified  experimentally. 

j 
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Appendix  A 

SIMPLIFICATION  OF  > 

The  coherency  of  a spherical  acoustic  wave  propagating  vertically 
into  the  lower  atmosphere  can  be  described  by  Its  first  and  second 
statistical  moments.  These  moments.  In  turn,  depend  on  the  perturbation 
[||'(^)  - 'I'q(?‘)]  0^  Its  complex  phase.  Consequently,  It  can  be  observed 
from  Eq.  (3.28)  that,  to  compute  these  moments.  It  Is  necessary  to 
compute  ^ )({;*(?"  )^  . Therefore,  (?' (r")^  Is  calculated 

in  this  appendix  for  spherical  wave  propagation.  From  Eqs.  (3.30)  and 
(3.31) 


'*'l(r)=^  “ ' h ^-rs-s-r-  d^'r 


V 


1 


(A.l) 


B^(r',r"),  the  covariance  function  of  the  acoustic  refractive  index,  is 


B„(f',r")  = <n^(f')n,(r")> 


Then 


/ . V -ik  (r'-r")  /’/*B  (?,,?„ 

<^l(r')^l(r")>  ® JJ  TV 

Vs  12 

exp[ikg  I (ri-r2)+|r'-r^  |-|^"■^2l  } ^ 

|r'-?^  Mr"-?2l 

• d^r^d^r2 


(A.2) 


(A. 3) 
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where  = (xp  * (x^.  P*'^**’*^  approx Iniat ion 

[Strohbehn,  1968],  and  forward  scatter; 


\rT-2 

= \/x  +y^  = 


X + 


li 


kM  -ik, 
x'«"e  > 


,.2  y.2 

•^a  g 

Jx"^  ’ ?jF. 


H. 


Bn(Fpr2) 


''s  ''s 


x^X2(x'-x^l{x"-X2)®’'^ 


ik. 


2 2 
^gl  y.2 
2^7  ■ IT^ 


ly^-yg/  l^u  *yg2l^ 
* TTx'-yT" 


d^r^d^r^  (A. 4) 


The  representation  of  r as  a quadratic  function  of  y^^  is  valid 
only  for  small  y^.  But  for  larger  values  of  y^  the  integrand 
oscillates  rapidly.  Hence  the  integration  over  y^  coordinates  has 
been  extended  to  infinity.  The  two  dimensional  spectral  representation 
of  B^(fpr2)  in  a locally  homogeneous,  isotropic  atmosphere,  with 
smooth  variations  in  vertical  direction  as  in  Eqs.  (2.10e),  (2.13)  and 
(2.16)  is 


B^(ri,r2) 


0 


e d < 


u 


(A. 5) 


Substituting  the  above  expansion  in  the  expression  for  (r' (r")^. 


. 2 

Ik  x,  X 

(x’-xj  -Ik  -y'  J- 

2k  x*  ' “ “ 

d 


2 

K^Xr 


'?0^vx  -x2;-iK^-y^  3J, 
a 


(A.9) 


(A.IO) 
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-Ik  •{?  ^ - y"  d^K  dx  dx, 
a ra  X ■'a  X > 02' 

. 


(A.ll) 


X ; x'.x"  = 6x  ; -^.n  ; x^-x^  » C 


yap  = ^ 5 ^am^K'ya  • Kl  ♦ Ul  5 Lg. 


'am  •'a  'a  ’ n’ 


the  outer  scale  of  turbulence.  6x  is  of  the  order  of  few  A ; its 

O 


maximum  value  being  the  number  of  cycles  in  the  acoustic  pulse  times 


A,.  Contributions  to  the  expression  in  the  above  equation  come  from 

a 


lower  values  of  also  |6x/x|  « 1.  Hence  terms  that  include  C 


and  6x/x  in  the  exponent  will  be  neglected.  Extending  the  limits  of 


integration  on  C to  « 


X 00  00 


exp  -<5|- «Mx)'  ■ (A-'2) 


Since 


00 

yF„(n,5. 


(A. 13) 


A r 

= 2»kj  J \) 


exp 

a 


- n .2-  . 

IK  .y  — d K dn 

a ■'am  x a 


(A. 14) 


Since  the  turbulence  is  isotropic,  and 


i 

J 
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(A.15) 


JL 


cv 


where  Jq(')  is  a Bessel  function  of  first  kind  and  zeroth  order. 


X » c 2 

<♦,(?' ■ 4,2kf  J Mpf-i^ex  (2)  I 

0 0 La  J 

• J^Ik  y 

01  a''am  x/  a 


(A. 16) 


For  a Kolmogorov  spectrum,  the  exponent  part  of  the  integrand  is  almost 
unity  over  the  useful  range  of  integrating  variables  for  typical 
values  of  acoustic  pulse  length.  As  a result 


X 00 

= 4.\2 1 ?)v'’ 

n n \ ' 


This  equation  is  used  in  computing  the  coherence  length  and  coherence 
function  of  vertical  acoustic  wave  propagation.  Evaluation  of 
<ki(x.  y^)|^>  is  as  follows.  From  Eq.  (A. 17), 


X oo 


(A. 18) 


Using  Eqs.  (2.20)  and  (2.23), 


X 00 


<k,(x,  y„)l^>  - (2i.k^)^J  J 0.033  C^Jf)  ' (A.t9) 

Oi-o'm 


a 11/6 


<l*,(x.  y„)|^>  = (0.0365)c2„k2A5'Vij  X 


(A. 20) 
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Appendix  B 

CONTRIBUTION  OF  DOUBLE-SCATTERING  TERM  TO  THE  MEAN  GREEN'S  FUNCTION 
OF  A SPHERICAL  ACOUSTIC  WAVE  PROPAGATING  IN  A TURBULENT  ATMOSPHERE 

To  compute  the  first  statistical  moment  of  the  acoustic  wave.  It 
is  necessary  to  determine  the  contribution  of  the  double-scattering 
term  in  Eq.  (4.17)  to  the  mean  of  a spherical  acoustic  wave  propagating 
in  a turbulent  atmosphere.  The  second  term  in  the  diagranmatic  series 
in  Eq.  (4.17)  is  defined  as  G,„(r,0)  I(x,y  ),  and  I(x,y  ) is  com- 

oO  d 

puted  in  this  appendix.  The  mathematics  is  the  same  as  in  Appendix  A, 
and  like  algebra  follows  in  Eq.  (4.18). 


2 2 


= (nr)  yjC  )",('-2)> 

X2>x^ 

ikadr-r^l  + |r2-r^  I + r^  - r) 
|r-r2l  1^2-^! I 

Using  the  paraxial  approximation 


(B.l) 


Kx.V 


= X 


ik)  ff 


B„(r,.r2) 


Xj>X|  (x-XjXxj-x,)*, 


• exp 


ik. 


ir  TT  |2 


2(x^o<2^  2(x2-x^  ) * 2x^ 


■'g 

2x 


d^r^d^2 


(B.2) 


From  Eqs.  (2.10e),  (2.13)  and  (2.16) 
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Note  that  depends  on  Ixpx^l.  From  Eqs.  (B.2)  and  (B.3) 


3) 


I(x.y^)  = X 


(-»-) 

'2n'  JJ  JJ  Tx-x^yCx^-x^llTj 


0 0-“ 


• I3(x.y^)d2..^dx,dx2 


where 


a I 2(x-x,)  ^(x^-x^ ) 2x^  2x  ) 


- i<  • (y  i-y  o) 
a al  oi2' 


Because 


ik  lll!a2l?ali!  + -7 

^a  I 2(x2-x^)  2x^  j ^,,1 


2- 


d y 


al 


'2TTix,  (xp-x, ) ' n 

L 2 a L 


’'^a^a2 


2x, 


2k, X, 

a c 


- K^a2 


then. 


■2Tiix,  {x,-x, ) 1 r y^  ik'^x,  (x,-x,  )1 

■ L ’ «P  S - 

, t d J L a c 


2 . 2 


1%^  * .2(-'  * “'y,: 
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(B.4) 


(B.S) 


(B.6) 


(B.7) 


APPENDIX  C 

SIMPLIFICATION  OF  I(x,  0/2,  y^) 


To  determine  the  intensity  of  the  RASS  received  signal,  it  is 
necessary  to  compute  which  is  the  ensemble  average  of  Ie^T 

A, 

This  quantity  is  derived  by  the  simplification  of  I(x,  0/2,  y^)  in 
Eq.  (5.38). 


'<«•»«•  ' / / "f--  >'am'  «:■» 


D'  and  D"  are  circles  of  radii  x'0/2  and  x"0/2,  respectively. 
Repeating  for  convenience,  the  definitions  given  in  Chapter  V, 


X - (X'  , x")/2;  = (y;  - y;;  );  )/2  (C.2) 


d'  = x'  e/2;  d"  = x"  6/2;  and  d ■ x 6/2  (C.3) 

a a a 


«(y„)  ■ 


1 y s d 
■^a  a 

0 > *^3 


(C.4a) 

(C.4b) 


the  above  integral  can  be  written  as 


l(x.  e/2,  y^„)  = f f ll(y;)  W(y;)  H(x,y^^)  d^  y^  d^  j;  (C.5) 


Making  the  change  of  variables  to  y^^^  and  y^p  results  in 


fflacimiw  fioi  MX  wumed 


This  double  integral,  is  the  overlap  area  of  two  circles  of  radii 
d'  and  d"  . Since  x'  = x";  A (x,§/2,  y ) is  evaluated  as  in 

a a d Olin 

Fried  [1967] 


1/2  4d  2 cos'^  !y  /(2d  )!  -y  (4d^-y^ 

I a / •^am  a | -^atn'  a -^atm 

y •-  2d  (C.8) 
■^am  a 

n y > 2d 

0 •'am  a 

The  maximum  value  of  the  above  expression  is  vrd^  . ^a(^’®/2,y^i)/(Trdg) 

is  plotted  in  Fig.  13.  From  Eg.  (C.6),  I (x,9/2,y^jjj)  does 


y /(2d  ) 
'am  0 


FIG.  13.  THE  OVERLAP  INTEGRAL 
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APPENDIX  D. 

COMPUTATION  OF  POWER-REFLECTION  COEFFICIENT 
IN  THE  PRESENCE  OF  A LINEAR  TEMPERATURE  PROFILE 

The  mathematics  involving  the  computation  of  received  power  in  an 
atmosphere  with  a linear  temperature  profile  is  considered  in  this 
Appendix.  It  can  be  observed  from  Eq.  (5.49)  that  the  phase  of  the 
acoustic  wave  has  a quadratic  dependence,  and,  therefore,  use  is  made 
of  Fresnel  sine  and  cosine  integrals.  The  quantity  F [defined  in 
Fqs.  (5.49)  and  (5.50)]  is  evaluated 


? -i(2k  r-a) 

F = I cos[k^r  - nr  ] e ^ dr  (D.l) 


n = agk^/(4T^)  . (D.Z) 

and  a is  a constant.  Initially,  the  Fresnel  sine  and  cosine  integrals 
are  defined  as  follows  [Abramowitz  and  Stegun,  1968], 


X 


Cf(x) 

= 1 cos  (J  t^) 

dt 

(D.3a) 

0 

X 

S^(x) 

= / (? 

dt 

(D.3b) 

0 


Some  useful  properties  of  these  integrals  are 


Using  these  properties,  F is  evaluated  to  be 


where 


F = J- 


cos(n  +01  )[C.(A  ) - C,(B  )] 
m '■  f m f m 


* Sin(n  R‘  * a,[S^(A^)  . S,(BJ] 


+ cos(n  Rp  -a)[C^(A  ) - C^(B  )] 


+ sin(n  Rp  - a)[S^(A  ) - S^(B  )] 


+ i I cos(n  R^  +a)[-S^(A^)  + S^(B^)] 


+ sin(n  R +oi)[Ca(A  ) - C (B  )] 
' m f'  m'  f'  m'-" 


+ cos(n  R^  - a)[S^(Ap)  - S^(Bp)] 


sin(ri  - a)[C^(Ap)  - C^(Bp)] 


(D.7) 


e = 


a A /T 
g a 0 


2X  / 

2A  \ 

2^.  / 

2A  . 

R = -T-  f 

1 § ) ; R = 

— ( 

1 + — ^ 

m 3 ' 

Ap  ^ P 

6 V 

? rp 

= -n  (Ri  - R„);  B„  = -r, 

m ^ v I m'  m ^ -rr 


(R?  - Rm) 

TT  2 m 


V ■ Oi- 


(D.8) 

(D.9) 

(D.IO) 

(D.ll) 
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For  ag  = 10  °K/km  and  300°K;  [R^l  > |(2A3)/6|  = |(2TQ)/ag|  = 

60  km.  Hence  |Rp|  is  much  larger  than  R^  and  R2  for  ranges 


the  acoustic-wave  train  is  short,  and  R,  - R,  = n X , then 

I c dC  a 


T\^X^/2  which  is  the  maximum  value  of  |F|.  This  maximum  value 

oC  a 


of  |F|  occurs  when  2X  /X  = [1  - 3x/(2A  )].  The  modified  Bragg- 

d 6 o 


scatter  condition  for  the  coherent  reflection  of  electromagnetic  energy 


from  an  acoustic  pulse  at  x for  |6x/(2X^)[  « 1,  therefore,  is 


1 + ^ 
0 


(D.23) 


which  can  be  explained  by  the  fact  that  peak  power  is  reflected  back 
when  the  interaction  between  electromagnetic  and  acoustic  waves  is 
maximum. 


1 


i- 

i 
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APPENDIX  E 

COMPUTATION  OF  NORMALIZED  RECEIVED  POWER 
IN  THE  PRESENCE  OF  HORIZONTAL  WINDS 


The  mathematical  details  involving  the  computation  of  received 

power  in  the  presence  of  constant  horizontal  winds  are  presented  in  this 

appendix.  An  asymptotic  expression  for  the  received  power  is  derived, 

using  an  asymptotic  expression  for  a Bessel  function  of  the  first  kind 

and  zeroth  order  J-(x).  Equations  (5.66)  and  (5.67)  define  F as 
0 w 

a quantity  proportional  to  the  strength  of  the  received  signal 


( cos  S^(r')  e 


-i(2k  r'-a)  , 

® d-^r’ 


(E.l) 


S„(f)  = k/ 


I ? 

L S J 


where 


r = (cosO,  cos  r,  sine,  sin  ^ sine) 

<»>=  (<"«><«,>.<“;>) 

<Wx^=  0;  = constant;  = constant 

<“y>'  ”h  *«  • <“2>  ' “h 
and 

then  from  Eq.  (E.2) 

r W 

= •'a''  ^ 


(E.2) 


(E.3a) 

(E.3b) 


(E.3c) 

(E.3d) 


(E.4) 


Some  useful  definitions  and  properties  of  the  Bessel  function  of  first 
kind  and  zeroth  order  are  listed  below  [Watson,  1966]: 
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Let  X be  the  height  of  the  acoustic  pulse.  Because  |sine|  < 1, 

k X w,/c  <<  1 implies  that  the  effect  of  horizontal  winds  can  be 
a h 0 ^ 

neglected  and  F = F as  in  Eq.(  5.68).  For  k xew,/c^  >>  1 
^ w wmax  ^ ' a h 0 

and  0 < 1,  a large-argument  approximation  for  Jq(')  in  Eq.  (E.9) 
can  be  used.  This  approximation,  however,  is  not  good  over  the  entire 
range  of  0 because  the  argument  of  the  Bessel  function  becomes  0 for 
0 = 0.  This  approximation  gives  rise  to  a minor  error  in  the  final 
result  because  the  integrand  has  sin  0 as  a factor.  In  addition, 
sin  6 fa  for  sin  (0/2)  is  small;  therefore 


Using  Erdelyi’s  theorem  [Popoulis,  1968J,  I-{r')  is  evaluated 

9 

asymptotically  to  be 


Note  that  the  sinusoidal  term  in  this  equation  varies  little  over  the 
acoustic  pulse  when  n^(,(w^/c^)(e  is  much  less  than  unity.  Under 
the  Bragg-scatter  condition  = 2X^  and  using  Eqs.  (5.68),  (E.IO)  , 
and  (E . 12) , therefore 


^ cos^Ek^x  0 w^/(2Cq)  - 3V4] 


wmax 


(k^x  e 


Because 


P /P 
rw' 


rmax 


|fj'/F 


2 

wmax 


(E.13) 


(E.14) 


where  is  the  power  received  in  presence  of  horizontal  winds  and 

^rmax’  9iven  by  Eq.  (5.41),  is  the  power  received  in  a static  atmosphere. 
Therefore 


where 


P^w  ^ 20.37  cos^Ek^x  6 w^/(2c^)  - 37r/4] 
^rmax  (k^x  0 w^/c^^)^ 

6 Vs  ""  ■'  • 


(E.15) 
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Naval  Weapons  Center 

Attn:  Code  5515,  H.  F.  Blazek 

China  Ixikc , CA  93555 

Dr.  H.  L,  Blood 
Technical  Director 
Naval  Undersea  Center 
San  Diego,  CA  95152 

Naval  Research  Laboratory 
Attn:  Code  5200,  A.  Brodzinsky 
1555  Overlook  Avenue,  SW 
Washington,  D.C.  20375 

Naval  Research  Laboratory 
Attn;  Code  7701,  J.  D.  Brown 
1555  Overlook  Avenue,  SW 
.Vashington,  D.C.  20375 


Naval  Research  Laboratory 
Attn:  Code  2627,  Mrs.  D.  Folen 
4555  Ovei’look  Avenue,  SW 
Washington,  D.C.  20375 

Dr.  Robert  R.  B'ossum 
Dean  of  Research 
Naval  Postgraduate  School 
Monterey,  CA  93940 

Dr.  G.  G.  Gould 
Technical  Director 
Naval  Coastal  System  Laboratory 
Panama  City,  FL  32401 

Naval  Ocean  Systems  Center 
Attn:  Code  7203,  V'.  E.  Hildebrand 
271  Catalina  Boulevard 
San  Diego,  CA  92152 

Naval  Ocean  Systems  Center 
Attn:  Code  753,  P.  M.  Johnson 
271  Catalina  Boulevard 
San  Diego,  CA  92152 
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Donald  E.  Kirk 
Professor  and  Chairman 

Electronic  Engineer,  SP-304 
Naval  Postgraduate  School 
Monterey,  CA  93940 


Naval  Research  laboratory 
Underwater  Sound  Reference  Division 
Technical  Library 
P.  0.  Box  8337 
Orlando,  FL  32806 


Naval  Air  Development  Center 
Attn:  Code  01,  Dr.  R.  K.  Lobb 
Johnsville 

Warminster,  PA  18974 

Naval  Research  laboratory 
Attn:  Code  5270,  B.  D.  McCombe 
4555  Overlook  Avenue,  SW 
Washington,  D.C.  20375 

Capt . R . B . Meeks 
Naval  Sea  Systems  Command 
NC  *3 

2531  Jefferson  Davis  Highway 
Arlington,  VA  20362 

Dr.  H.  J.  Mueller 
Naval  Air  Systems  Command 
Code  310 
JP  J»1 

1411  Jefferson  Davis  Highway 
Arlington,  VA  20360 

Dr.  J.  H.  Mills,  Jr. 

Naval  Surface  Weapons  Center 
Electronics  Systems  Department 
Code  DF 

Dahlgren,  VA  22448 

Naval  Ocean  Systems  Center 
Attn:  Code  702,  H.  T.  Mortimer 
271  Catalina  Boulevard 
San  Diego,  CA  92152 

Naval  Air  Development  Center 
Attn:  Technical  Library 
Johnsville 

Warminster,  PA  18974 

Naval  Ocean  Systems  Center 
Attn:  Technical  Library 
271  Catalina  Boulevard 
San  Diego,  CA  92152 


Naval  Surface  Weapons  Center 
Attn:  Technical  Library 
Code  DX-21 
Dahlgren,  VA  22448 

Naval  Surface  Weapons  Center 
Attn:  Technical  Library 
Building  1-330,  Code  WX-40 
White  Oak 

Silver  Spring , MD  20910 

Naval  Training  Equipment  Center 
Attn:  Technical  Library 
Orlando,  FL  32813 

Naval  Undersea  Center 
Attn:  Technical  Library 
San  Diego,  CA  92152 

Naval  Underwater  Systems  Center 
Attn:  Technical  Library 
Newport , RI  02840 

Office  of  Naval  Itesearch 
Electronic  and  Solid  State 
Sciences  Program  (Code  427) 
800  North  Quincy  Street 
Arlington,  VA  22217 

Office  of  Naval  Research 
Mathematics  Program  (Code  432) 
800  North  Quincy  Street 
Arlington,  VA  22217 

Office  of  Naval  Research 
Naval  Systems  Division 
Code  220/221 
800  North  Quincy  Street 
Arlington,  VA  22217 

Director 

Office  of  Naval  Research 
New  York  Area  Office 
715  Broadway,  5th  Floor 
New  York,  NY  10003 
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Office  of  Naval  Research 
San  Francisco  Area  Office 
One  Hallidie  Plaza,  Suite  601 
San  Francisco,  CA  94102 

Director 

Office  of  Naval  Research 
Branch  Office 
495  Summer  Street 
Boston,  MA  02210 

Director 

Office  of  Naval  Research 
Branch  Office 
536  South  Clark  Street 
Chicago,  IL  60605 

Director 

Office  of  Naval  Research 
Branch  Office 
1030  East  Green  Street 
F^sadena , CA  91101 

Mr.  H.  R.  Riedl 

Naval  Surface  Weapons  Center 

Code  WR-34 

White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Naval  Air  Development  Center 
Attn:  Code  202,  T.  J.  Shoppie 
Johnsville 

Warminster,  PA  18974 

Naval  Research  Laboratory 
Attn:  Code  5403,  J.  E.  Shore 
4555  Overlook  Avenue,  SW 
Washington,  D.C.  20375 

A.  L.  Slafkovsky 
Scientific  Advisor 
Headquarters  Marine  Corps 
MC-RD-1 

Arlington  Annex 
Washington,  D.C.  20380 

Harris  B.  Stone 

Office  of  Research,  Development, 
Test  and  Evaluation 
NOP-987 

The  Pentagon,  Room  5D760 
Washington,  D.C.  20350 


Mr.  L.  Sumney 

Naval  Electronics  Systems  Command 
Code  3042,  NC  #1 
2511  Jefferson  Davis  Highway 
Arlington,  VA  20360 

David  W.  Taylor 
Naval  Ship  Research  and 
Development  Center 
Code  522.1 
Bethesda , MD  20084 

Naval  Research  Laboratory 
Attn:  Code  4105,  Dr.  S.  Teitler 
4555  Overlook  Avenue,  SW 
Washington,  D.C.  20375 

Lt . Cdr.  John  Turner 
NAVMAT  0343 
CP  «5,  Room  1044 
2211  Jefferson  Davis  Highway 
Arlington,  VA  20360 

Naval  Ocean  Systems  Center 
Attn ; Code  746 , H . H . Wieder 
271  Catalina  Boulevard 
San  Diego,  CA  92152 

Dr.  W.  A.  Von  Winkle 
Associate  Technical  Director 
for  Technology 

Naval  Underwater  Systems  Center 
New  London,  CT  06320 

Dr.  Gernot  M.  R.  Winkler 
Director,  Time  Service 
US  Naval  Observatory 
Massachusetts  Avenue  at 
34th  Street,  NW 
Washington,  D.C.  20390 


Other  Government  Agencies 

Dr . Howard  W . Etzel 

Deputy  Director 

Division  of  Materials  Research 

National  Science  Foundation 

1800  G Street 

Washington,  D.C.  20550 
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Mr.  J.  C.  French 
National  Bureau  of  Standards 
Electronics  Technology  Division 
Washington,  D.C.  20234 

Dr.  Jay  Harris 
Program  Director 
Devices  and  Waves  Program 
National  Science  Foundation 
1800  G Street 
Washington,  D.C.  20550 

Los  Alamos  Scientific  laboratory 
Attn:  Reports  Library 
P.  0.  Box  1663 
Los  Alamos,  NM  87544 

Dr . Dean  Mitchell 
Program  Director 
Solid-State  Physics 
Division  of  Materials  Research 
National  Science  Foundation 
1800  G Street 
Washington,  D.C.  20550 

Mr.  F.  C.  Schwenk,  RD-T 
National  Aeronautics  and 
Space  Administration 
Washington,  D.C.  20546 

M.  Zane  Thornton 
Deputy  Director,  Institute  for 
Computer  Sciences  and  Technology 
National  Bureau  of  Standards 
Washington,  D.C.  20234 


Nongovernment  Agencies 
Director 

Columbia  Radiation  I^iboratory 
Columbia  University 
538  West  12oth  Street 
New  York,  NY  10027 

Director 

Coordinated  Science  laboratory 
University  of  Illinois 
Urbana , IL  61801 


Director  of  laboratories 
Division  of  Engineering  and 
Applied  Physics 
Harvard  University 
Pierce  Hall 
Cambridge,  MA  02138 

Director 

Electronics  Research  Center 
Hie  University  of  Texas 
Engineering-Science  Bldg.  112 
Austin , TX  78712 

Director 

Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 

Director 

Electronics  Sciences  Laboratory 
University  of  Southern  California 
Los  Angeles , CA  90007 

Director 

Microwave  Research  Institute 
Polytechnic  Institute  of  New  York 
333  Jay  Street 
Brooklyn,  NY  11201 

Director 

Research  laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Cambridge , MA  02139 

Director 

Stanford  Electronics  laboratory 
Stanford  University 
Stanford , CA  94305 

Stanford  Glnzton  laboratory 
Stanford  University 
Stanford , CA  94305 

Officer  in  Charge 
Carderock  laboratory 
Code  18  - G.  H.  Gleissner 
David  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda , MD  20084 
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Dr.  Roy  F.  Potter 
3868  Talbot  Street 
San  Diego,  CA  92106 


